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FOREWORD 


This  report  was  prepared  under  Task  1\  0130901,  IR-IS9.  It  Is  concerned 
with  the  complicated  relations  and  Interdependence  of  the  various  critical 
parameters  as  well  as  relations  between  the  results  of  various  shock 
sensitivity  tests.  The  same  material  was  presented  in  a  series  of  eight 
lectures  In  1980-81.  It  should  be  of  interest  In  the  areas  of  explosive 


sensitivity  and  propellant  safety. 


jL~,  9/2/ — 

JAMES  F.  PROCTOR 
By  direction 
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CRITICAL  PARAMETERS  FOR  DETONATION  PROPAGATION  AND  INITIATION  OF  SOLID  EXPLOSIVES 

INTRODUCTION 

It  is  the  purpose  of  this  report  to  examine  the  available  relevant 
Information  and  to  elucidate  the  very  complicated  relationships  and  Inter¬ 
dependence  of  the  various  critical  parameters  for  detonation  propagation  and 
Initiation  In  solid  HE.  In  particular,  the  role  of  test  diameter  relative  to 
critical  diameter  In  affecting  the  critical  pressure  for  initiating  detonation 
(and  vice  versa)  will  be  described  and  Illustrated.  Comparisons  will  also 
be  made  between  experimental  values  of  shock  sensitivity  from  different  gap 
tests  and  the  wedge  test.  Critical  density  (in  the  form  of  %  TMD),  critical 
particle  size,  and  critical  temperature  for  propagation,  and  critical  pressure 
for  shock  initiation  will  also  be  considered. 

Since  there  will  necessarily  be  cross  references  to  critical  diameter  and 
gap  test  values,  these  will  be  defined  here,  and  later  described  and  discussed 
In  separate  sections.  The  critical  diameter  (dj  is  defined  for  a  cylindrical 
solid  charge  as  the  minimum  diameter  at  which  a  steady  state  detonation  can 
propagate.  The  gap  test  Is  a  widely  used  measure  of  shock  sensitivity.  In  Its 
simplest  form,  it  consists  of  a  cylindrical  donor,  attenuator,  and  acceptor. 

We  shall  consider  only  those  tests  in  which  the  shock  attenuator  (gap)  is  a 
solid  material.  The  test  arrangement  is  that  of  Fig.  1.  In  a  series  of  test 
firings,  the  gap  thickness  is  varied  until  the  50$  point  Is  reached,  l.e.  that 
gap  thickness  *or  which  detonation  occurs  In  50%  of  the  trials.  That  50$ 
point  gap  is  defined  as  the  test  value.  If  the  test  has  been  calibrated,  the 
test  value  can  also  be  given  as  shock  pressure  at  the  end  of  the  gap  (Pg). 

If  the  Hugoniot  of  the  explosive  acceptor  is  known,  this  can  be  converted  to 
Pc1  the  transmitted  pressure  which  is  the  minimum  Initiating  or  critical  pressure 
to  Initiate  detonation  In  that  acceptor  In  that  test. 

To  achieve  steady-state  detonation,  a  cylindrical  charge  must  have  d  >  dc 
and  the  initiating  shock  must  have  an  amplitude  P  >  P  Together  these 
conditions  are  necessary  and  sufficient  for  steady-state  detonation.  This 
report  will  consider  each  in  turn. 
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CRITICAL  DIAMETER 

Historical 

Sometimes  the  critical  diameter  is  defined  as  the  smallest  diameter  of  a 
charge  at  which  propagation  of  an  undamped  detonation  wave  is  still  possible. 
This  definition  has  been  avoided  here  because  it  includes  the  case  of  gases 
and  liquids  in  which  a  constant  velocity  front  appears,  but  the  propagation  is 
by  means  of  fading  and  re-ignition.  In  other  words,  the  steady-state  detonation 
structure  is  absent.  At  the  present  time,  this  phenomenon  has  not  been  observed 
in  solids  and  it  is  believed  that  the  sustained  front  in  detonating  solid  HE  has 
the  Zeldovitch-von  Neumann-Doering  (ZND)  model  structure.  According  to  this 
model,  the  detonation  wave  with  one-dimensional  flow  has  the  structure: 
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The  plane  of  the  von  Neumann  shock  Is  separated  from  the  Chapman-Jouguet  (C-J) 
plane  by  a  distance  called  the  reaction  zone  length.  Reaction  can  start  at  the 
shock  front;  it  must  be  completed  at  the  C-0  plane.  The  pressure  drop  between 
these  two  planes  results  from  a  combination  of  heterogeneous  (hot  spot)  and 
homogeneous  reactions;  In  the  latter,  the  von  Neumann  pressure  Is  maintained 
for  tha  induction  time  (i.e.  throughout  most  of  the  reaction  zone)  after  which 
the  pressure  falls  very  rapidly  to  the  C-J  value.  In  the  case  of  heterogeneous 
reaction,  the  pressure  starts  to  drop  immediately  after  the  von  Neumann  shock 
and  reaches  the  C-J  value  with  a  less  steep  gradient. 

The  structural  unit  bounded  by  the  C-J  and  von  Neumann  planes  is  propagated 
unchanged  in  a  steady  state  detonation.  On  the  other  hand,  the  gas  products 
follow  the  C-J  plane  in  a  continuously  changing  structure  as  the  wave  propagates. 

^cording  to  the  definition  chosen,  the  critical  initiating  pressure  for 
cv' indr'  al  charges  at  the  critical  diameter  is  then  the  von  Neumann  spike 
ssure,  and  the  critical  initiating  pressure  pulse  is  that  between  the 
jn  Neuiranr  and  C-J  planes  of  the  detonation  front.  These  two  planes  are 
separated  by  3  distance  a  (reaction  zone  length)  and  a  time  (reaction  time). 

/..other  .ay  of  saying  this  is  that  at  the  critical  diameter,  a  cylindrical 
charge  can  be  cu*  in  two  cylinders  of  half  length.  One  half  can  be  used  as 
tab  :lcnor,  the  other  as  an  acceptor  in  a  gap  test.  The  50?  point  should  then 
be  zc  0  gap. 

Fj.'  essentially  voidless  materials  of  95%  or  greater  theoretical  maximum 
density  (TMD),  .:e  detonation  velocity  at  critical  conditions  (Dc)  is,  of  course, 
Igs  han  the  Ideal  detonation  velocity  (D^).  However,  it  is  not  a  great  deal 

lecs  at  high  %  T.iD.  For  example,  pressed  and  cast  TNT  at  98%  TMD  have  D_/Dj 

1  L  1 
values  of  0.95  and  0.94,  respectively  .  DNT  at  99%  TMD  has  a  D./D,  value 

of  0.85".  On  the  other  hand,  high  porosity,  granular  charges  exhibit  much 

Vrice,  D. ,  "Shock  Sensitivity,  A  Property  of  Many  Aspects,"  5th  Symposium 
(International)  on  Detonation,  ONR  ACR-184,  U.S.  Gov.  Printing  Office, 
Wash'ngton,  DC,  1972;  pp  207-217. 

2 

Price,  D. ,  Erkman,  J.  0.,  Clairmont,  Jr.,  A.  R. ,  and  Edwards,  D.  J.,  "Explosive 
Characterization  of  Dinitrotoluene,"  Combust.  Flame,  Vol .  14,  1970,  pp  145-148. 
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lower  values;  e.g.,  TNT  at  50%  TMD  has  a  D_/D .  value  of  0.65*  and  25y  AP  at 

3  ci 

51%  TMD,  about  0.45  .  These  values  are  given  for  illustration  o>:ly,  not  for 
generalization.  D  ,  like  d  ,  is  a  function  not  only  of  the  chemical  composition 
of  the  charge,  but  also  of  its  initial  particle  size  (6),  temperature,  and  %  TMD. 
In  the  conventional  1-D,  plane  theory,  the  C-J  pressure  ( Pj )  is  given  by 

Pj  a  p0D?/(y+1)  (1) 

This  theory  does  not,  of  course,  consider  diameter  effect;  that  can  be  treated 
only  with  a  2-  or  3-D  model  with  a  reaction  zone  of  finite  thickness.  However, 
we  can  use  Eq.  1  to  obtain  estimates  of  the  orders  of  magnitude  involved. 

Thus  at  Dc/D^  «  0.95,  (Pj)c~  0*9  Pj-  In  the  case  of  high  density  TNT, 

PvN  =  1.25  Pj4.  Hence  (PvN)c~  1.25(Pj  )c  ~  0.9  x  1.25  ?.  »  1.125  x  1894  ■ 

213  kbar. 

5  6 

Khariton  ’  evidently  first  attempted  to  express  the  diameter  effect  on 
detonation  velocity  D  and  the  critical  diameter  d  quantitatively.  Like  earlier 

v 

workers,  he  explained  the  effect  by  energy  losses  caused  by  rarefactions  entering 
the  reaction  zone  from  the  lateral  surfaces.  And  the  principle  he  proposed  was: 

dc~  2tc  (2) 

where  c  is  the  average  sound  velocity  in  the  detonation  products  and  x  is  the 
detonation  reaction  time.  Thus  xc~  distance  traveled  (one  radius)  by  the 
rarefaction  before  it  quenches  the  reaction. 

15 - 

Price,  D. ,  Clairmont,  Jr.,  A.  R. ,  and  Jaffe,  I.,  "Explosive  Behavior  of 
Ammonium  Perchlorate,"  Combust.  Flame,  Vol.  11  No.  5,  1967,  pp  415-425. 

Also  ibid.  Vol.  13,  No.  1,  1969,  pp  104-108. 

4 

Coleburn,  N.  L.  and  Llddiard,  Jr.,  T,  P.,  "Hugoniot  Equations  of  State  of  Several 
Unreacted  Explosives,"  J.  Chem.  Phys.,  Vol.  44,  No.  5,  1966,  pp.  1929-1936. 

5Khariton,  Yu.  B.  and  Rozing,  V.  0.,  Dokl  .  AN  SSSR,  Vol.  26,  1940,  p.  360. 

Through  Ref.  7  and  8. 

^Khariton,  Yu.  B.,  Collection:  "Problems  of  the  Theory  of  Explosives,"  Izd-vo 
AN  SSSR,  Moscow,  1947.  Through  Ref.  7  and  8. 
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From  Eq.  2,  Khar  1  tor;  derived 


for  the  same  explosive  under  different  conditions,  e.g.  different  Initial 
particle  size.  He  also  derived 


(3) 


(4) 


i!is  derivation  is  described  in  Ref.  9,  but  an  estimate  can  be  obtained  from 
the  1-D  theory.  There 


c  =  (y/y+1)D  (5) 

Hence  c  D  and  Eq.  2  produces  Eq.  4.  The  latter  must  be  applied  cautiously, 

however,  because  Y  and  hence  the  proportionality  function  can  differ  from 

explosive  to  explosive;  Y  is  also  a  function  of  the  density. 

Khariton's  principle  stimulated  a  great  deal  of  experimental  work  on  dr 

7  c 

In  Russia.  But  since  the  advent  of  the  electromagnetic  method  to  measure  mass 
velocity  as  a  function  of  time,  and  to  estimate  x,  it  has  been  recognized  that 
Eqs.  2-4  are  certainly  not  quantitative^’10  and  cannot  be  used  to  estimate 
even  relative  values  of  t. 

^Dremin,  A.  N. ,  Savrov,  S.  D. ,  Trofimov,  V.  S. ,  and  Shvedov,  K.  K. ,  "Detonation 
Waves  in  Condensed  Media,"  Izd-vo  Nauka,  Moscow,  1970.  Translation  AD  751417. 

Q 

Dremin,  A.  N.  and  Trofimov,  V.  S. ,  "Nature  oi  the  Critical  Detonation  Diameter 
of  Condensed  Explosives,"  Combustion,  Explosion,  and  Shock  Waves,  Vol .  5, 

No.  3,  1969,  pp.  208-12. 

g 

Rybanin,  S.  S. ,  "Effect  of  Inhomogentles  in  an  Explosive  on  the  Critical 
Detonation  Diameter,"  PMTF,  >'ol .  10,  No.  4,  1969,  pp  141-144,  through  a 
translation  journal. 

10Kravtsov,  V.  V.  and  Silvestrov,  V.  V. ,  "Effect  of  Low  Temperature  on  Detonation 
Parameters  of  Cast  Trotyl,"  Combustion,  Expl .  and  Shock  Waves,  Vol.  15,  No.  3, 
1979,  pp.  387-390. 
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11  12 

Slightly  later,  Jones  In  England  and  Eyrlng  et  al.  In  USA  Introduced 
analytical  expressions  for  D  as  a  function  of  d.  The  latter,  for  unconfined 
charges.  Is 


(6) 


where  "a"  Is  the  reaction  zone  length,  assumed  a  constant.  This  Is  a  linear 

relationship  between  the  ratio  of  detonation  velocities  end  d’1.  It  fits 

D(d)  data  of  most  explosives  well  at  large  d  and  not  at  all  as  d  approaches 
13 

dc  .  Jones'  relationship  is  more  complicated;  It  fits  the  data  somewhat 
better,  but  not  well 13. 

14 

Recently  Campbell  and  Engel ke  have  added  a  correction  term  to  Eq.  6 
to  obtain: 


A 

m 


(7) 


where  A  and  E  are  length  parameters  and  R  Is  the  charge  radius.  Eq.  7  fits 
the  D(d)  data  of  more  than  a  dozen  essentially  voidless  explosives  down  to 
d  *  d  . 


Jones,  H.,  "The  Theory  of  the  Dependence  of  the  Rate  of  Detonation  of  Solid 
Explosives  on  the  Diameter  of  the  Charqe,"  Proc.  Roy.  Soc. ,  Vol.  A189,  1947, 
p.  415. 

12Eyring,  H. ,  Powell,  R.  E.,  Duffey,  G.  H. ,  and  Parlln,  R.  H.,  "The  Stability 
of  Detonation,'1  Chem  Revs.,  Vol.  45,  1949,  p.  69. 

13Mal1r»,  K.  E.  and  Campbell,  A.  W. ,  "°art1cle-S1ze  Effects  in  Explosives  at 
Finite  and  Infinite  Diameters,"  J.  Appl.  Phys.  Vol.  28,  1957,  pp.  63-69. 

Also  Mai  in,  Campbell  and  Mautz,  "Partlc'le  Size  Effect  in  One-  and  Two-Component 
Explosives,"  2nd  ONR  Symposium  on  Detonation,  1955,  pp  360-381. 

^Campbell,  A.  W.  and  Engelke,  R.,  "The  Diameter  Effect  in  High  Density 
Heterogeneous  Explosives,"  Sixth  Symposium  (International)  on  Detonation, 

ACR  221,  ONR,  U.S.  Gov.  Printing  Office,  Washington,  1978;  pp.  642-663. 
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Dremln,  et  al.7,8  considered  that  Jones'  nozzle  theory11,  Eyrlng's  curved 
front  theory12,  and  subsequent  work1^’19  based  on  "extremely  simplified  schemes 
of  flows"  have  produced  little  or  no  advance.  In  other  words,  a  rigorous  theory 
of  detonation  limits  for  condensed  explosives  has  not  yet  been  developed. 

O 

Dremln  formulated  the  problem  and  described  the  fundamental  physical  principles 

which  he  believed  could  be  the  basis  for  constructing  such  a  theory.  His 

discussion  is  the  most  helpful  reference  for  failure  limits  encountered  in  the 

literature.  Another  hopeful  initial  treatment  has  been  recently  published  by 

20 

Chalken  and  Edwards  . 

Measurement 

The  critical  diameter  can  be  measured  In  several  ways.  The  explosive  can 

be  cast  or  machined  into  a  right  cylindrical  cone.  When  It  is  initiated  at  the 

base,  d.  Is  the  cone  diameter  at  which  detonation  fails.  The  detonation  can 
C  21 

be  followed  optically  or  electronically  .  This  method  has  the  disadvantage 
that  propagation  from  larger  to  smaller  diameters  along  the  cone  results  In 

overboosting  which  may  show  quenching  at  a  diameter  somewhat  less  than  the 
tf - 

Wood,  W.  W.  and  Kirkwood,  J.  G.,  "Diameter  Effect  in  Condensed  Explosives. 

The  Relation  between  Velocity  and  Radius  of  Curvature  of  the  Detonation  Wave," 
J.  Chem.  Phys.,  Vol  22,  No.  11,  1954,  pp  1920-23. 

16Cook,  M.  A.,  "The  Science  of  High  Explosives ," Relnhold,  New  York,  1958. 

17Dubnov,  L.  V.,  "Losses  in  a  Detonation  Wave,"  Russian  J.  Chem.  Phys.,  Vol.  30, 
No.  10,  1960, pp  1124-5. 

18Evans,  M.  W. ,  "Detonation  Sensitivity  and  Failure  Pianeler  in  Homogeneous 
Condensed  Materials,'  J.  Chem.  Phys.,  Vol.  36,  No.  1,  1962,  p.  193. 

19Rempel ,  G.  G. ,  collection:  "Explosives"  No.  52/9,  39,  through  Ref.  7. 

20Cha1ken,  R.  F.  and  Edwards,  John  C.,  "A  Kinetic  Lattice  Approach  to  Detonation 
In  Heterogeneous  Explosives,"  Sixth  Symposium  (International)  on  Detonation, 
ACR  221,  ONR,  U.S.  Gov.  Print.  Office,  Washington,  D.  C. ,  1978;  pp  344-351. 

21Jaffe,  I.  and  Price,  D. ,  "Determination  of  the  Critical  Diameter  of  Explosive 
Materials,"  ARS  Journal,  Vol  32,  1962,  pp  1060-1065. 
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true  dc.  Another  method  Is  to  use  a  stepped  cylinder  with  long  enough  segments 
to  allow  any  overboosting  to  die  out.  At  d  »  d  ,  a  charge  length  i  >  3d  Is 

v 

quite  adequate;  but  as  d  ■+  dc,  the  Russians  use  i  of  10-12  d.  Finally,  a 
series  of  cylinders  of  varying  diameter  can  be  used,  and  D(d)  determined.  Then 
Eq.  7  can  be  used  to  select  d  In  the  Interval  between  the  diameter  at  which 
detonation  failed  (dc-)  and  the  next  larger  d  where  It  propagated  (dc  ).  This 
last  method  seems  the  best  although  the  value  of  dc  so  obtained  may  differ  very 
little  from  the  average  of  dc+  and  dc‘  If  the  experimental  increment  (dc+  -  dc") 

Is  small. 

There  are  Insufficient  data  in  the  literature  to  compare  conical  and 
cylindrical  determinations  of  d  from  the  sane  charge  by  the  same  investigators. 

v 

If  the  angle  of  the  cone  Is  kept  very  small.  It  Is  doubtful  that  the  difference 
In  the  method  of  measurement  could  be  distinguished  from  the  effect  of  differences 
In  charge  preparation  that  would  Inevitably  occur  in  two  different  laboratories. 
Consequently,  no  distinction  has  been  made  In  the  method  used  to  obtain  the  d 

w 

values  listed  In  the  later  tables. 

Because  it  Is  difficult  to  prepare  charges  as  cylinders  c'  very  small 
diameter,  It  Is  much  simpler  to  measure  a  critical  layer  height  of  an  explosive 
slab.  For  a  confined  layer  this  thickness  was  designated  hc  and  its  relation 
to  d£  was  given  as 

d  -  2hc  (8) 

??  1 
by  Belyaev  and  Sukoyan  for  a  flat  plate  of  thickness  h  of  TNT  (pQ  -  1.57  g/cm‘ , 

6q  *  400  -  800y);  the  explosive  was  confined  on  the  underside  by  a  5  ran  thick 

witness  plate  of  lead.  It  Is  reasonable  to  exoect  some  such  relationship 

between  d.  and  h  .  but  how  general  and  how  qualitative  It  may  be  Is  an  open 
c  c  23 

question.  In  a  recent  LASL  report  ,  the  results  of  critical  wedge  thickness 
tests  are  given;  these  tests  too  had  a  confining  witness  plate  (brass)  as  well 
as  some  lateral  confinement.  Data  for  four  cast  and  four  pressed  explosives 

72 - 

Belyaev,  A.  F.  and  Sukoyan,  M.  k. ,  "Detonability  of  Some  Explosives  with 
Increase  In  External  Pressure,"  Comb.,  Exp!.,  and  Shock  Waves,  Vol.  3,  No.  1, 
1967,  pp  11-13. 

23 

Urizar,  M.  J.,  Peterson,  S.  W. ,  and  Smith,  L.  C.,  "Detonation  Sensitivity 
Tests,"  LA-7193-MS,  April  1978. 
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on  which  d  values  are  also  available  are  given  In  Table  1.  Since  the  ratio 
c 

d  / h  varies  from  1.2  to  7.6,  It  does  not  appear  that  Eq.  8  or  a  similar  one 
c  c 

can  now  be  used  for  anything  except  establishing  an  order  of  magnitude.  However, 
a  more  controlled  stujy,  e.g.  Identical  preparation  and  charges  (pQ,  composition, 
particle  size)  might  Improve  this  situation. 

In  addition  to  being  a  function  of  the  chemical  composition  of  the  charge, 
the  value  of  the  critical  diameter  Is  a  strong  function  cf  1)  degree  of 
compaction/porosity,  2)  Initial  particle  size  of  a  granular  charge  or  the  nature 
of  the  casting  of  non-granular  charge,  and  3)  the  test  temperature.  These 
variables  will  be  considered  In  turn. 

Poroslty/%  TMD  Effect  (Critical  Density) 

If  we  represent  the  theoretical  maximum  density  by  TMC,  then  the  percent 
porosity  Is  (100  -  %  TMD).  For  TNT-Ilke  materials,  dc  decreases  as  %  TMD 
increases  (see  Fig.  2);  for  materials  like  ammonium  perchlorate  and  Its  mixtures, 
dc  Increases  a:>  %  TMD  Increases,  (See  Fig.  3).  These  curves,  dividing  the  d-% 

TMD  plane  into  detonation  and  failure  areas,  are  called  detonablllty  curves. 

For  any  given  HE,  at  a  specific  %  TMD,  there  Is  a  critical  diameter  for  detonation, 
d  .  Similarly  for  any  fixed  charge  diameter,  there  Is  a  critical  value  of  %  TMD 

V 

(or  pQC)  for  propagation  of  detonation.  Hence  one  may  speak  of  either  critical 

diameter  or  critical  density.  Other  critical  parameters  described  later  must 

also  be  considered  In  parallel  although  these  two  are  probably  the  most  studied. 

At  the  time  data  of  Figs.  2  and  3  were  determined,  their  trends  were  used 

27 

to  classify  HE  Into  Group  1  (TNT-Ilke)  and  Group  2  (AP-lIke)  explosives  . 

However,  further  work  revealed  that  there  was  ari  upturn  In  the  Group  1  curve 
at  high  %  TMD,  and  In  the  Group  2  curve  at  low  %  TMD.  In  other  words,  the 
general  failure  curve  appears  to  be  U  shaped  although  both  branches  are  not 
generally  found  In  easily  accessible  experimental  ranges.  Fig.  4  shows  the 
reversals  for  NQ  (Group  1)  at  80-85%  TMD  and  AP  (Group  2)  at  40-25%  TMD. 

Materials  showing  Group  2  behavior  at  high  %  TMD  are  those  likely  to  dead-press, 
l.e.  become  subcrltlcal  at  high  density. 

^Prlce,  D. ,  "Contrasting  Patterns  in  the  Behavior  of  High  Explosives," 

Eleventh  Symposium  (International)  on  Combustion,  Combustion  Inst., 

Pittsburgh,  1967;  pp  693-702. 
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TABLE  1 

Critical  Heights  and  Diameters  for  Various  Charges 


isiupni'i 


a  lues 
d_ 


Material 

(g/cm3) 

hc  (mm) 

(mm) 

c 

Ref. 

•  | 

Cast  Charges 

\ 

i 

50/50  Pentollte  c 

1.700 

1.39 

6.7 

21 

4.8  | 

Comp  B(A)  c 

1.713 

1.42 

4.3 

14 

3.0  | 

Comp  B-3  c 

1.728 

0.82 

6.2 

21 

7.6 

75/25  Cyclotol  c 

1.752 

1.51 

6.0 

14* 

4.0  ; 

Pressed  Charges 


0ATB 

1.724 

0.73 

5.3 

25c 

’  7.3  1 

P8X  9404 

1.78-1.85 

0.46 

1.18 

14 

2.6  j 

Comp  A- 3 

1.63 

0.56 

<2.2 

24 

<3.9 

TNT 

1.63 

?.16d 

2.62 

14 

i.*  ! 

! 

*Ref.  23. 

If  density  range  was 

short,  values 

of  replicate 

runs 

1 

were  averaged;  ] 

large  ranges  in  pQ  are  Indicated. 

a.  Value  for  77/23  cyclotol, 

a 

b.  Value  for  po  ■  1.8  g/crrT  . 
d.  Pressed  at  65°C. 


DiAMETER  d  (mm) 
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DETONATOR 


sm 

DONOR 

m?PJ 

ACCEPTOR 

FIGURE  1  SIMPLE  GAP  TEST 


FIGURE  2  DETONABILITY  LIMITS  IN  THE  CHARGE  DIAMETER- 
POROSITY  PLANE,  GROUP  1.  EXAMPLE  SHOWN  IS  TNT 
(GRAIN  SIZE  0.07  TO  0.2  mm)  DATA  (REF.  26) 
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PERCENT  THEORETICAL  MAXIMUM  DENSITY  (%  TMD) 


F'G'JRE  3  DETONABILITY  LIMITS  IN  THE  CHARGE  DIAMETER 
POROSITY  PLANE,  GROUP  2.  EXAMPLE  SHOWN  IS  AP 
(AVERAGE  PARTICLE  SIZE  10  M)  DATA  (REF.  3).  WHERE 
VALUE  HAS  BEEN  BRACKETTED,  (+)  INDICATES  JETONAT!ON, 
(-),  FAILURE 


MINIMA  IN  DETONABILITY  CURVES.  <  •  DETONATION;  O  FAILURE).  (REF.  28) 
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The  compaction  effect  on  d  results  In  rather  different  D(d)  vs  %  TMD 

patterns  for  the  two  HE  groups  as  shown  in  Fig.  5.  These  patterns  as  well 

as  the  detonability  curves  (dc  vs  %  TMD)  can  be  used  for  classification. 

Although  the  extrapolated  portions  of  the  D  vs  %  TMD  curves  cf  Fin.  5a  (Group  1) 

are  drawn  to  show  a  marked  curvature  as  they  approach  the  failure  threshold, 

such  curvature  may  be  difficult  to  detect  experimentally  for  the  rapidly 

reacting  charges  of  Group  1.  Thus  TNT  appears  to  have  linear  curves  even  to 

29 

the  failure  limit  except  for  the  smallest  diameter  charge  tested  ;  its 
detonability  plot  showed  definite  curvature. 

Figs.  5  a  &  b  illustrate  very  clearly  that  Group  1  materials  approach 
Ideal  behavlcr  (D^  and  small  reaction  zones)  at  high  %  TMD  whereas  Group  2 
HE  become  more  ideal  at  low  %  TMD.  The  different  behaviors  are  attributed 
to  the  dominance  of  homogeneous  reaction  in  Group  1  materials  and  of  hetero¬ 
geneous  reaction  In  Group  2  materials  as  discussed  in  Ref.  27. 

Relatively  few  detonability  curves  are  available,  and  not  a  great  many 
Investigations  of  D(d)  have  been  carried  out.  From  the  available  data,  members 
of  Group  1  are: 

TNT 

All  the  more  powerful  common  organic  explosives* 

Mixtures  of  the  above  with  each  other  and  inerts 
NQ** 

DNT*** 


29 

Stesik,  L.  N.  and  Akimova,  L.  N. ,  "An  Indirect  Method  of  Estimating  Reaction 
Zone  Width  of  a  Detonation  Wave,"  Russ.  J.  Phys.  Chem. ,  Vol .  33,  No.  8,  1959, 
pp  148-151. 

*So  classified  because  no  anomalies  In  D  measurements  have  been  reported.  If 
one  should  occur,  the  specific  HE  should  be  investigated  in  detail. 

**Group  2  behavior  at  high  densities  (Figs.  4  and  8). 

***Group  2  behavior  at  high  densities  (Fig.  8).  Ref.  27  quoted  a  Russian 
statement  that  "weak  Individual  explosives  e.g.  dinitro* compounds  of  aromatic 
hydrocarbons"  exhibited  Group  2  behavior.  Subsequent  study  showed  this  false 
for  DNT 2  pressed  from  material  of  s*10y. 
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AP 

AN 

Their  mixtures  with  fuel 
HN31 

32 

Ammonium  Azide 
Hydrazonium  Azide33 

Probably  Lead  Azide  in  view  of  above  and  fact  that  it  dead  presses. 
Particle  Size  (6)  Effect  (Critical  Particle  Size) 

The  particle  size  effect  on  dc  is  the  same  for  both  groups  of  explosives: 
the  smaller  6,  the  smaller  d  .  Hence  decreasing  6  increases  the  ability  of  the 
charge  to  propagate  detonation.  Fig.  C  Illustrates  this  for  TNT;  Fig.  7,  for 
AP.  rig.  8  is  also  Included  to  show  Group  1  and  Group  2  plots  on  the  same  scale. 

The  failure  curve  of  10  -  80u  AP  In  Fig.  7  comes  from  reference  34.  The  AP 
used  In  that  work  was  obtained  from  a  sieve  cut;  It  had  a  microscopically 
observed  particle  size  range  of  10  -  80y,  but  the  particle  size  distribution 
and  median  were  not  determined.  The  weight  median  particle  size  would  of  course 
be  smaller  than  the  mean  of  the  projected  area  diameter  of  microscopic  observations. 
Hence  it  Is  not  surprising  that  this  curve  lies  between  those  for  the  weight 

OC 

median  particle  sizes  of  10  and  25y.  More  recent  work  gives  a  detonability 
curve  for  5y  AP  which  lies  below  all  the  others  in  Fig.  7. 


3*Pr1ce,  D. ,  Liddiard,  T.  P.,  Jr.,  and  Drosd,  R.  D. ,  "The  Detonation  Behavior 
of  Hydrazine  Mononitrate,"  NOL  TR  66-31,  Apr  1966. 

32Yakovleva,  G.  S.,  Kurbangalina,  P.  Kh . ,  and  Stesik,  L.  N. ,  "Detonation 
Properties  of  Ammonium  Azide,"  Combust.,  Exp!.,  and  Sh.  Waves,  Vol.  13, 

No.  3,  1977,  pp  405-7. 

33Yakovleva,  G.  S. ,  Kurbangalina,  R.  Kh.,  and  Stesik,  L.  N.»  "Rate  of  Detonation 
of  Hydrazonium  Azide,  ibid,  Vol.  10,  No.  2,  1974,  pp  270-4. 

^Gor'kov,  V.  A.  and  Kurbangalina,  R.  Kh.,  "Some  Data  Concerning  the  Detonation 
Ability  of  Ammonium  Perchlorate,"  Combust,,  Explosions,  and  Shock  Waves, 

Vol.  2,  Mo.  2,  1966,  p.  12. 

35Akimova,  L.  N.  and  Stesik,  L.  N.,  "Detonation  Capacity  of  Perchlorate  Explosives," 
Combust.  Expl.  and  Shock  Waves,  Vol.  12,  No.  2,  1976,  pp  247-251. 
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FIGURE  6  DETONATION  FAILURE  LIMIT  CURVES  FOR  TNT.  (REF.  26  &  28! 
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The  only  available  data  I  have  encountered  for  dr  vs  5  at  a  fixed  p  (or  %  TMD) 
are  those  of  Gor'kov  and  Kurbangalina  .  Fig.  9  snows  their  data  for  56.4% 

TMD  AP  (pQ  *  1.1  g/cm^).  Our  data  for  10  and  25v  AP  fall  very  close  to  the 
same  curve.  As  drawn,  the  curves  appear  to  level  off  with  dc~110  ircn  at 
6  s,  500y.  This  va1ue  of  500y  may  be  inexact;  It  is  possible  that  there  Is  a  low 
gradient  it  continuing  particle  size  effect  on  very  coarse  AP.  The  existence 
of  curves  such  as  Figs.  6-9  implies  that  for  a  given  Initial  temperature  and 
density,  there  Is  a  critical  &  above  which  steady  state  detonation  cannot 
propagate. 

The  examples  of  AP  charges  were  at  low  %  TMD.  For  charges  of  granular 
TNT  or  other  HE  easily  detonable  at  high  %  TMD,  Fig.  6  shows  that  the  Initial 
particle  size  effect  is  almost  negligible  at  high  %  TMD.  This  would  be  expected 
since  pressing  to  voidless  density  tends  to  eliminate  differences  In  Initial  <5. 
Melt-Casting  Effect 

~  IP 

Cybulsk!  et  al.  gave  an  excellent  example  of  the  effect  of  different 
casting  techniques  (and  consequent  different  grain-sized  products)  on  the 
critical  diameter  of  TNT.  Their  results  were: 


'asting  Method 


Poured  clear 
Creamed* 


dc,cm 

3.17  <  dc 
2.20  <  dc  <  2.54 
1.26  <  dc  <  1.66 


Creamed*  plus  10%  fines  1.26  <  d  <  1.66 

V 

Examples  given  in  Ref.  1  showed  that  the  more  rapid  the  cooling  rate  of  the 
melt,  the  smaller  the  grain  size  and  the  lower  the  d  value. 

V 

Temperature  Effects  (Critical  Temperature) 

Temperature  would  be  expected  to  affect  the  reaction  rate  and  hence  the 
critical  diameter.  Such  an  expected  temperature  gradient  is  shown  In  Fig.  10 
for  gr  ^ular  TNT  .  It  is  for  a  charge  of  ofl  *  1.0  g/cm"  [and  presumably 
fixed  fi]  over  a  temperature  range  of  -l^i0  to  75°  C,  and  the  limit  curve 

-jg  ■ 

Cybulski,  W.  B  Payman,  W.,  and  Woodhead,  D.  W. ,  "Explosion  Waves  and 
Shock  Waves  VI*.  The  Velocity  of  Detonation  in  Cast  TNT,"  Froc.  Roy.  Soc. 
(London),  Vol.  1S7A,  1949,  pp  51-72. 

37 

Belyaev,  A.  F.  and  Kurbangalina,  R.  Kh.,  "Effect  of  Initial  Temperature  on 
Critical  Diameter  of  Nitrogylcerin  and  Trinitrotoluene,"  Russian  J.  Phys. 
Chem. .  Vol,  34,  No.  3,  1960,  pp  285-9. 

*Term  as  used  at  NSWC  i.e.,  melt  stirred  until  it  becomes  milky  because  of 
formation  of  small  crystals;  at  this  point  it  is  poured  into  the  mold. 


NSWC  TR  80-339 


divides  the  plane  Into  super-  and  sub-critical  regions.  Thus  for  a  fixed 
Initial  density  and  particle  size,  there  exists  a  critical  temperature  (t  ) 

77  v 

below  which  detonation  cannot  occur. 

Although  t  has  not  been  featured  In  the  literature  to  the  extent  of  d  , 

w  W 

It  Is  a  very  Important  critical  parameter  and  most  evident  In  the  explosive 
behavior  of  propellants  and  propellant  components.  Thus  AN  prills  at  pour 
density  are  subcrltlcal  under  heavy  confinement  (ID  of  38.1  mm)  until  their 
temperature  Is  Increased  to  140°C  or  higher38.  AP  also  shows  a  marked  effect 
of  temperature  with  a  gradient  of  about  0.1  mm/°C^  as  compared  to  the  TNT  gradient 
of  about  0.02  mrr/°C  (Fig.  10).  Several  double  base  propellants,  tested  In  the 
confinement  of  the  NOL  large  scale  gap  test,  failed  to  detonate  at  -60°C  and 
-32°C  but  did  detonate  at  25°C  and  66°C38. 

In  addition  to  the  usual  temperature  effect  on  dc,  described  above,  there 
exists  the  possibility  of  a  temperature  effect  on  the  structure  of  cast  materials; 
this  physical  change  can  then  cause  a  change  in  dr<  Cast  TNT,  when  rapidly 
cooled  to  77.4°K  and  allowed  to  return  to  its  Initial  temperature  of  290°K 
by  slow  heating,  showed  no  visible  change  in  Its  microstructure.  Nevertheless, 

d  had  been  decreased  from  15  to  7  mm  (fine  crystal  casting)  and  from.  27.5  to 

^  3 

9  mm  (coarse  crystal  casting)*;  the  initial  p.  was  reduced  by  0.01  g/cm  from 

3  10  0 

1.62  and  1.60  g/cm  ,  respectively  .  Although  two  such  temperature  cycles 
equal lized  the  d  values  of  the  two  castings  at  6.5  mm,  further  cycling  had  no 
additional  effect.  In  other  words,  after  the  temperature  cycling,  the  cast 
material  had  a  dc  value  still  about  2.5  times  that  of  pressed  TNT  at  the  same 
density. 

Chemical  Composition  Effects 

In  the  earlier  Russian  literature,  It  was  reported  that  a  mixture  of  a 

Group  1  HE  with  a  Group  2  HE  (50/50  TNT/AN)  exhibited  a  maximum  In  Its  d  vs 

35  ^ 

%  TMD  curve.  Subsequent  Russian  work  makes  this  seem  doubtful.  Fig.  11 

shows  detonablllty  curves  for  a  series  of  Group  1/Group  ?  mixtures  -  In  this  case, 

HMX/AP.  Each  composition  gave  a  smooth,  monotonic  detonabillty  (d  vs  pQ)  curve 

38 

Price,  D. ,  Jaffe,  I.,  and  Roberson,  G.  E.,  "Shock  Sensitivity  of  Solid 
Explosives  and  Propellants,"  Ind.  chlm.  Beige,  1967,  No.  32  (Spec.  No.), 
pp  506-510. 

*The  dc  value  measured  at  77.4°K  was  9  mm  In  both  cases. 
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with  Group  2  behavior  predominating  at  90%  or  more  AP,  10%  or  less  HMX;  Group  1 

behavior,  at  85%  or  less  AP,  15%  or  more  HMX.  We  found  a  similar  shift  In  the 

detonablllty  curves  from  two  mixtures;  I.  62.5/18.75/18.75  AP/Al/Wax  and 

II.  50/15/15/20  AP/Al/Wax/HMX  or  80/20  I/HMX.  Fig.  12  shows  the  shift  from 

Group  2  behavior  to  that  very  close  to  Group  1  caused  by  the  addition  of  20% 

HMX.  Whereas  15%  HMX  sufficed  to  convert  an  AP/HMX  mixture  to  Group  1  behavior 

35 

(Fig.  11),  about  32%  TNT  Is  required  for  such  a  conversion  .  However,  the 
change  from  one  behavior  to  the  other  is  a  continucus,  smooth  variation  following 
the  variation  in  composition. 

If  the  two  component  mixture  is  made  up  of  members  of  the  same  group,  the 
mixture  will  also  exhibit  a  detonablllty  curve  typical  of  that  group.  The 
magnitude  of  its  d£  will  depend,  of  course,  on  its  composition.  Fig.  13 
shows  a  dc  vs  composition  curve  for  cyclotols  (RDX/TNT),  with  linear  variation 
of  dc  between  the  value  for  RDX  at  pQ  *  1.0  p/cm3  to  TNT  at  pQ  ■  1.0  g/cm3. 

Figure  12  showed  that  addition  of  A1  to  AP  depressed  the  curve  ( 1 . e . 
decreased  dc)  but  did  not  change  Its  form.  This  may  also  be  the  case  for 
porous  TNT  charges^.  Dlmza's  data^  were  obtained  at  a  constant  volume 
fraction  of  HE  and  no  value  for  the  unadulterated  HE  was  given;  hence,  It  Is 
impossible  to  separate  the  effect  of  the  additive  from  the  effect  of  the 
change  in  porosity  It  caused.  It  should  be  noted  here,  however,  that  addition 
of  A1  to  cast  TNT  definitely  decreases  d  from  2.7  cm  to  0.7  cm  at  20%  Al^*. 

v 

Addition  of  water  to  AP  or  TNT  Increases  d  ;  the  effect  Is  very  much 
14  ^ 

greater  for  AP  ,  possibly  because  of  its  solubility.  Water  behaves  In  both 

cases  as  an  inert  diluent.  Addition  of  wax  to  TNT,  RDX,  or  HMX  Increases  dc 

{In  this  case,  dlmza's  data  show  a  monotonic  Increase  or  no  change  with 

simultaneous  increase  in  diluent  and  decrease  In  porosity.  Since  the  latter 

change  decreases  dc  of  Group  1  HE,  it  Is  evident  that  the  wax  must  have  acted 

as  an  inert  diluent.)  On  the  other  hand,  addition  of  wax  to  AP  decreases  dg; 

see  Fig.  7.  This  Is  typical  of  any  AP/fuel  mlxture3^’3^.  Both  this  fact  and 

35 

the  fact  that  the  decrease  Is  much  greater  if  the  fuel  is  finely  divided 

™D1mza,  G.  V. ,  "Critical  Detonation  Diameter  for  an  Explosive  Containing  an  Inert- 
Additive,"  Comb.,  Expl.  and  Shock  Waves,  Vol.  12,  No.  2,  1976,  pp  244-47. 
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are  Indicative  of  an  oxidant-reductant  reaction  contributing  to  the  detonation 
behavior.  A  further  indication  of  this  Is  the  Increase  in  D  observed  with 
the  addition  of  wax  to  A P3®. 

Finally,  It  should  be  noted  that  the  addition  of  inert  materials  to  voidless 

explosives  can  affect  d_.  (A1  in  cast  TNT  may  be  an  example  of  this.)  Apin 

41  c 

and  Stesik  added  relatively  dense  powdered  diluents  (CaCO^,  MgO,  BlgO^,  PbP, 
HgO,  and  W)  to  plastlsol  type  propellants.  They  found  that  such  diluents 
reduced  dc  according  to  the  formula: 


(n  +  nQ) 


T7T 


*  B 


*2 
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Price,  0.,  and  Clalrmont,  Jr.,  A.  R.,  "Explosive  Behavior  of  Simplified 
Propellant  Models,"  Combust.  Flame,  Vol.  29,  1977,  pp  87-93. 


(9) 


where  A  &  B  are  functions  of  the  diluent, 

n  =  concentration  cf  additive  particles, 

nQ  a  concentration  of  reaction  foci  in  matriy  before  addition  of  inert. 
The  first  term  of  Equation  9  is  the  average  distance  between  foci.  It  v/as 
assumed  that  the  Inert  particles  of  P>cmatr-jx  afforded  a  focus  where  shock 
reflection  and  intensification  could  produce  a  hot  spot. 

Very  few  data  have  been  available  for  the  effect  of  inert  diluents  In 
porous  explosives;  recently  such  data  have  been  compiled  and  supplemented 
Inert  diluents  in  HE,  at  concentrations  of  5%  or  greater  and  at  the  pour 
density  of  the  mixture,  either  show  no  effect  on  dc  (RDX/paraff in  up  to  ca. 
35%  paraffin)  or  increase  It  (TNT/talc,  paraffin,  or  NaCl  and  RDX/NaCl). 
Moreover,  curves  of  d  vs  s  (initial  particle  size  of  diluent)  showed  a 

v 

minimum.  In  other  words,  there  was  a  6  which  caused  the  smallest  Increase 
In  dc  at  a  fixed  concentration  of  inert. 


41 

Apin,  A.  Ya.  and  Stesik,  L.  N.,  "The  Chemical  Reaction  Mechanism  In  the 
Detonation  of  Explosives,  PMTF,  No.  2,  1965,  pp  142-145. 


^Shvedov,  K.K.,  Aniskin,  A.  I.,  II 'In,  A.  N.,  and  Dremin,  A.  N.,  "Detonation 
of  Highly  Diluted  Porous  Explosives  I.  Effect  of  Inert  Additives  on  Detonation 
Parameters,"  Comb.,  Expl .  and  Shock  Waves,  Vol.  16,  No.  3,  1980,  pp  324-331. 
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Confinement  Effect 

It  has  been  well  known  for  many  years  that  confining  a  charge  will  decrease 
dc*,  l.e.  the  charge  diameter  within  the  sheath  which  will  support  a  detonation. 
Very  few  quantitative  measurements  of  this  effect  have  been  made,  and  there  Is 
no  adequate  theoretical  treatment  with  which  predictions  can  be  made.  Apparently 
the  confinement  effect  Is  most  striking  for  materials  of  large  d  .  Thus  d  of 
AP  In  thin  cellophane  is  4-5  times  greater  than  that  in  glass  which  Is  In  turn, 

Oil 

4-5  times  greater  than  that  of  AP  In  copper  foil  .  Similarly,  ammonium  nitrate 
(AN)  at  pQ  ■  1.06  g/cnr  has  dc  i  100  mm,  i  80  mm,  and  <  3  mm  In  cellophane 
glass,  and  steel,  respectively^,  a  steel  wall  thickness  of  3  mm  was  adequate 
to  assure  detonation  In  a  4  mm  ID  steel  tube  filled  with  AN  at  pQ  *  0.8  g/cm3  ^3. 


*The  confinement  delays  the  expansion  of  the  gaseous  detonation  products  and  thus 
Increases  the  time  available  for  reaction  at  high  temperature  and  pressure  before 
a  rarefaction  reaches  the  reaction  zone.  In  this  manner,  the  charge  behaves 
like  a  much  larger  bare  charge.  Hence  confinement  has  increased  Its  "effective 
diameter." 

43 

Kurbangallna,  R.  Kh.  and  Patronova,  L.  I.,  "Effect  of  a  Steel  Sheath  on  the 
Critical  Detonation  Diameter  of  Condensed  Explosives,"  Comb, ,  Exp! . ,  and 
Shock  Waves,  Vol.  12,  No.  4,  1976,  pp  587-90. 
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CRITICAL  INITIATING  PRESSURE 


Background 

In  the  previous  section,  the  parameters  critical  for  propagation  of 
detonation  were  considered.  In  this  section  we  will  consider  the  parameter, 
critical  Initiating  pressure  (Pc1)  which  helps  to  define  the  threshold  for 
a  shock-to-detonatlon  transition  (SOT)  and  Is  an  inverse  measure  of  the  shock 
sensitivity  of  an  explosive. 

The  critical  Initiating  pressure  is  defined  as  the  minimum  pressure  (P^) 
required  to  Initiated  detonation  of  the  explosive  In  50%  of  the  trials.  It 
defines  a  threshold  for  the  Initiation  of  detonation  under  transient  conditions. 
The  value  of  P  ^  measured  In  any  given  test  Is  determined  not  only  by  tho 
chemical  composition  of  the  test  material  and  Its  physical  state  but  also,  In 
large  part,  by  the  physical  dimensions  of  the  test.  These  dimensions  control 
the  times  at  which  rear  and  lateral  rarefactions  can  arrive  at  the  shock  front 
and  relieve  the  pressure.  Thus,  each  test  devised  to  measure  Pc<  Is  apt  to 
use  a  different  pressure  pulse  (pressure-time  profile)  to  excite  reaction  of 
the  explosive.  Different  profiles  can  be  obtained  by  varying  the  composition 
and  shape  of  donor  explosives  In  the  gap,  booster,  and  wedge  tests  and  by 
varying  projectile  material,  shape,  and  velocity  In  the  projectile  impact 
and  wedge  experiments.  They  are  also  obtained  by  varying  the  diameter  of  the 
acceptor  explos4,,e.  The  explosive  response  differs  as  the  stimulus  differs, 
and,  hence,  gives  different  P  ^  values  In  different  tests.  Any  single 
adequately  designed  test  will  give  a  relative  sensitivity  rating  to  a  group 
of  different  HE  charges. 

44 

Some  time  ago  It  was  shown  theoretically  that  a  limiting  pressure-time 

(P-t)  curve  exists  for  a  50%  chance  of  build-up  to  detonation  in  homogeneous 
TT 

Hubbard,  H.  W.  and  Johnson,  M.  H.,  "Initiation  of  Detonations,"  J.  Appl.  Phys. 

Vol .  30,  No.  5,  1959,  pp  765-69 
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explosives;  the  limiting  curve  divides  the  Impulse  plane  Into  a  detonation 
region  and  a  failure  region.  Brown  and  Whitbread46  demonstrated  by  Impact  of 
cylindrical  projectiles  on  two  physically  heterogeneous  explosives  that  the 
P^  determined  must  have  a  minimum  duration  for  detonation  to  occur.  They 
also  showed  qualitatively  that,  at  p.  pc,  and  a  duration  shorter  than  that 
required  at  the  critical  pressure  level,  detonation  could  also  be  initiated. 

Since  then  there  has  been  rather  general  agreement  that  initiation  of  detonation 
must  be  the  result  of  the  pressure-time  history  of  the  Initiating  shock,  and 
that  a  critical  limit  curve  such  as  Figure  14  must  exist  In  the  P-t  plane 
for  each  explosive. 

In  many  cases  of  explosive  loading  (e.g.,  a  gap  test),  maximum  shock 
pressure  can  be  determined  much  more  readily  than  the  pressure-time  profile. 

Hence,  many  of  the  shock  sensitivity  measurements  have  been  restricted  to 
pressure  measurements  only.  However,  It  Is.  generally  known  that  the  Pc1  will 
depend  on  Its  duration.  Consequently,  for  each  test  In  which  the  pressure-time 
history  of  the  shock  loading  Is  unknown,  the  experimental  details  are  carefully 
specified  (e.g.,  dimensions  and  character  of  the  donor  explosive  In  the  gap  test.) 

Solid  heterogeneous  explosives  (cast  and  pressed)  exhibit  the  following 
shock  sensitivity  behavior46,47:  a  critical  Initiating  pressure  (Pc1 )  for 
detonation,  breakout  of  detonation  downstream  from  the  shocked  boundary,  and 
a  consistent  decrease  of  both  the  run  length  and  delay  tine  to  steady-state 
detonation  with  increased  amplitude  of  the  applied  shock.  The  run  length 
Is  the  distance  from  the  plane  of  shock  entry  Into  the  explosive  to  the  plane 
In  which  steady-state  detonation  is  first  established.  The  delay  time  Is  the 
total  time  from  the  moment  of  shock  entry  to  the  time  at  which  steady-state 
detonation  begins.  Figure  15  illustrates  this  behavior;  In  particular,  It 
It - 

Brown,  S.  M.  and  Whitbread,  F.  G.,  "The  Initiation  of  Detonation  by  Shock 
Waves  of  Known  Duration  and  Intensity,",  Les  Ondes  de  Detonation,  Edition 
du  Centre  National  de  la  Recherche  Sclentifique,  Paris,  1962;  pp.  69-80. 

46 

Jacobs,  S.  J.,  Llddlard,  T.  P.,  and  Drlmmer,  B.  E.,  "The  Shock-to-Detonatlon 
Transition  In  Solid  Explosives,"  Ninth  Symposium  (International)  on  Combustion, 
Academic  Press,  N.Y.C.,  1963;  pp  517-529. 

47Pr1ce,  D.  and  Petrone,  F.  J. ,  "Dotonation  Initiated  by  High-Pressure  Gas 
Loading  of  a  Solid  Explosive,"  J.  Appl.  Phys.,  Vol.  35,  No.  3,  1964,  pp  710-714. 
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FIGURE  14  VARIATION  OF  INITIATING  PRESSURE  WITH  ITS  DURATION 


RUN  LENGTH/DELAY  TIME  TO  STEADY-STATE  DETONATION 

FIGURE  15  VARIATION  OF  RUN  LENGTH/DELAY  TIME  WITH  AMPLITUDE 
OF  INITIATING  SHOCK 
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shows  how  the  initiating  pressure  P^  (for  a  specific  test  charge)  decreases 
to  the  critical  Initiating  pressure  P^  as  the  run  length  increases  to  a  large 
value.  The  approach  to  Pc^  is  asymptotic.  If  the  charge  is  cylindrical  and 
at  its  critical  diameter,  the  value  of  P^  should  be  approximately  the  C-J 
pressure  for  that  explosive  at  d  s  d  ,  the  pressure  needed  just  to  sustain 

w 

steady-state  detonation. 

The  P^-diameter  relationship  should  follow  a  curve  such  as  Fig.  16  where 
the  upper  point  on  the  threshold  curve  is  at  (Pj)c  for  d  *  dc  (and  consequently 
at  a  time  duration  of  the  detonation  reaction  time  t  at  d  »  d  ).  The  larger 
the  diameter,  the  longer  the  path  and  travel  time  of  the  lateral  rarefaction 
waves  from  the  charge  periphery  to  the  charge  axis;  it  is  the  arrival  of  these 
release  waves  that  quenches  the  reaction.  Therefore  the  larger  d,  the  longer 
Pj  and  Tj  are  maintained.  Hence  the  reaction  zone  length  and  reaction  time 
decrease  with  increasing  charge  diameter  to  their  ideal  values.  As  d  increases 
above  dc,  P  ^  should  drop  sharply  and  approach  asymptotically  its  infinite 
diameter  value  where  t  also  assumes  its  ideal  value.  (A  similar  sharp  change 
in  P  ^  would  be  expected  in  approaching  the  critical  values  of  other  parameters, 
e.g.  tr,  6C<)  Fig.  16  indicates  that  the  test  value  of  Pc-  for  any  given 
explosive  will  depend  on  d/dc ,  as  it  is  known  to  do,  and  Ref.  48  gives  an  experi¬ 
mental  curve  P  .  vs  d  for  NEMO*  powder  at  1.58  g/cc  for  d  >  d  .  Inasmuch  as 

L>  I  v 

the  critical  diameter  value  is  a  function  of  temperature,  initial  particle  size, 

density,  composition,  and  confinement,  and  P  .  is  a  function  of  d/d it  follows 

cl  c 

that  P  ^  is  also  a  function  of  t,  6,  pQ,  composition,  and  confinement.  However, 

the  effect  of  these  variables  on  P  .  is  not  alone  through  their  effect  on  dc 

because  the  same  change  in  6,  pQ  and  composition  can  change  Pci  and  dc  in 

opposite  directions  as  we  will  demonstrate  later. 

Pci  can  be  measured  by  gap,  booster,  projectile  impact  and  wedge  tests 

provided  the  test  is  calibrated  and  Hugon-ot  data  for  the  non-reacting  explosive 

are  available.  Only  the  first  and  last  of  these  tests  will  be  considered  in 

detail.  We  believe  all  the  test  results  are  related  and  will  later  illustrate 

relations  between  gap  tests  and  between  gap  and  wedge  tests. 

_ 

Afanasenkov,  A.  N.  and  Voskoboinikov,  I.  M. ,  "Shock-Ignition  Sensitivity  of 
Balllstite,"  Comb.  Expl.  Shock  Waves,  Vol .  9,  Mo.  2,  1973,  pp  331-332. 

*A  NC/NG  double-base  powder  containing  40%  NG. 
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NOL  Large  Scale  Gap  Test 

The  NOL  large  scale  gap  test  (LSGT)  was  used  for  much  of  the  data  quoted 
later,  and  will  therefore  be  described  in  some  detail.  Figure  17  shows  that 
this  LSGT  is  a  conventional  gap  test  with  a  5  cm  dia  x  5  cm  long,  1.56  g/cc 
50/50  pentolite  donor  and  an  attenuator  gap  of  polymethyl  methacrylate  (PMMA) 
or  its  equivalent.  The  acceptor  explosive  ccr  be  tested  unconfined  or  in  the 
standard  moderate  confinement  of  a  steel  sleeve  (3,65  cm  ID,  4.76  cm  OD).  A 
mild  steel  witness  plate  Is  used  in  both  cases,  and  the  gap  length  is  varied 
until  the  50$  value  is  found,  that  is,  the  attenuation  at  which  a  hole  is 
punched  in  the  witness  plate  in  50%  of  the  trials.  The  system  donor/gap  has 
been  calibrated  [Fig.  18^]to  give  shock  pressure  as  a  function  of  gap  "'ecgth. 
This  50%  pressure  (P  )  at  the  end  cf  the  gap  can  be  converted  to  the  initiating 
pressure  (P  . )  transmitted  to  the  explosive  by  the  use  of  the  Hugoniot  of  the 
gap  material49  and  that  of  the  explosive,  e.g.,  cast  TNT4. 

The  Hugoniot  adiabat  or  equation  of  state  of  the  attenuator  PMMA  (p^  = 

1.18  g/cm  )  Is  required  for  the  calibration  cf  the  test.  Usually  this  is 
known  in  the  form  of 

U  =  aQ  +  su  (10) 

where  U,  u  are  shock  and  particle  velocity,  respectively,  and  aQ  and  s  are 
constants  specific  to  the  material.  Actually  PMMA  is  visco-elastic  and  its 
exact  Hugoniot  is  more  complicated  than  the  linear  relation  of  Eq.  10  (See 
Ref.  49).  However  Eq.  10  is  a  good  approximation  and  quite  adequate  for 
this  illustration.  In  combination  with  the  hydrodynamic  relationship  derived 
from  the  conservation  of  momentum, 

P  =  pqUu  (11) 

the  non-linear  relationship  between  pressure  P  and  u  in  the  shocked  PMMA  can 
be  obtained.  This  is  the  customary  and  convenient  relationship  with  which 
to  work  because  boundary  conditions  for  reflection  and  transmission  of  shock 
waves  at  the  interface  between  two  different  materials  require  that  pressure 

P  and  particle  velocity  u  must  be  equalized  across  the  boundary. 

_  - 

Erkman,  J.  0.,  Edwards,  D.  J.,  Clairmont,  A.  R. ,  Jr.,  and  Price,  D. , 
"Calibration  of  the  NOL  Large  Scale  Gap  Test;  Hugoniot  Data  for  Polymethyl 
MethacrylaLe,"  NOLTR  73-15,  Apr  1973. 
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Measurement  of  shock  velocity  vs  gap  thickness  Lg  or  of  particle  velocity 

vs  Lg  In  the  shocked  PMMA  together  with  the  P-u  Hugonlot  provides  the  Pg-Lg 

calibration  of  Fig.  18.  Both  U  and  u  measurements  were  made  on  the  longitudinal 

axis  and  compared  to  a  computation  of  the  hydrodynamic  flow  In  the  donor/PMMA 

49 

system  to  obtain  the  calibration  . 

Pg  is  generally  used  to  compare  different  HE  because  non-reactlve  Hugonlot 

data  are  not  available  for  all  explosives.  Where  such  data  are  available  the 

measured  P_  can  be  converted  to  P.j ,  the  critical  initiating  pressure  in  the 
9  C1  4 

explosive.  For  example  the  cast  TNT  Hugoniot  and  that  for  PMMA  are  shown 

In  Fig.  19.  Pg  at  the  50%  point  was  44  kbar  in  PMMA  (pQ  *  1.18  g/cm^).  At 

point  X  (L_)  the  P,u  values  are  44  kbar  and  0.92  mm/us,  respectively.  At  the 
9  o 

interface  with  the  denser  TNT  (1.61  g/cm  ),  a  reflected  shock  (dashed  line)  is 

sent  back  Into  the  PMMA  to  raise  Its  pressure  to  52  kbar  which  Is  also  the 

pressure  of  the  shock  transmitted  to  the  TNT  (point  b).  Fig.  20  shows  a 

comparison  of  P  with  P  .  for  six  cast  and  three  pressed  HE  all  near  voidless 
60  9  c‘ 

density  .  There  Is  no  reversal  in  the  sensitivity  ordering  although  the  P^ 

are  16-30%  greater  than  the  corresponding  Pg.  Hence  the  Pg  rating  shrinks 

but  does  not  distort  the  true  sensitivity  scale  based  on  P 

The  standard  donor  of  the  LSGT  is  approximately  point  Initiated;  hence, 

the  detonation  front  In  the  donor  is  spherical  and  the  hydrodynamic  flow  In 

the  gap  and  In  the  acceptor  Is  divergent.  The  transmitted  pressure  falls  off 

rapidly  as  Fig.  18  shows.  From  optical  observations  of  the  tetryl*/PMMA  system, 

51 

the  shock  pressure  decreases  to  half  Its  Initial  amplitude  In  about  one  usee  . 
Walker  and  Wasley52  estimated  that  such  a  pulse  was  equivalent  to  a  square 
pulse  of  an  amplitude  of  0.9  P^  and  1.6  usee  duration.  Although  a  one-dimensional 
hydrodynamic  flow  computation  on  the  tetryl/PMMA  system  produced  a  half-width 

^Price,"  D. ,  "Large  Scale  Gap  Test:  Interpretation  of  Results  for  Propellants," 
NAVWEP  7401,  Mar  1961. 

5*Jaffe,  I.,  Toscano,  J.,  and  Price,  D. ,  "Behavior  of  Plexiglass  under  Shock 
Loading  by  a  Tetryl  Donor,"  NOLTR  64-66,  2  Jul  1964. 

52Walker,  F.  E.,  and  Wasley,  R.  J.,  "Critical  Energy  for  Shock  Initiation  of 
Heterogeneous  Explosives,"  Explosl vstoffe,  Vol.  17,  No.  1,  1969,  p.  9. 

3 

♦Pressed  tetryl  pellets  at  1.51  g/cm  were  originally  used  as  the  LSGT  donor 
in  place  of  the  present  standard  pentolite. 
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GAP  THICKNESS,  Lg  (IN.x  100) 


FIGURE  18  CALIBRATION  CURVE  FOR  NOL  LSGT  WITH  PENTOLITE  DONOR 


P  (kbar) 
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FIGURE  19  DETERMINATION  OF  Pci  FROM  Pg  IN  CAST  TNT 
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PRESSURE  AT  END  OF  GAP,  Pg  (kbar) 

FIGURE  20  COMPARISON  OF  SHOCK  LOADING  AT  50%  POINT  WITH  CRITICAL 
INITIATING  PRESSURE  IN  CHARGE.  (REF.  50) 
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of  4  vs  ,  analogous  two-dimensional  computations  Indicated  a  value  of 

about  2  vs,  In  fair  agreement  with  the  original  estimate  from  experimental 

work. 

Like  all  such  tests,  the  LSGT  exhibits  a  diameter  effect.  Inasmuch  as 

doubling  the  length  of  the  standard  booster  Increases  L  but  does  not  affect 

54  9 

the  measured  pressure  P  ,  rear  rarefactions  affect  the  duration  of  the 
pressure  pulse  produced  by  the  standard  donor,  but,  from  the  standpoint  of 
Initiation,  the  pulse  Is  already  of  Infinite  duration.  It  Is  therefore 
the  lateral  rather  than  the  rear  rarefaction  that  has  the  more  Important 
effect  In  reducing  the  pressure  amplitude.  Because  confinement  of  the  charge 
delays  the  arrival  of  lateral  rarefactions  at  the  charge  axis,  we  find,  for 
the  same  measured  P^,  that  the  diameter  of  the  equivalent  unconfined  charge 
Is  approximately  twice  the  core  diameter  of  the  confined  charge,  l.e.,  76  mm. 

As  In  the  case  of  dc,  both  Pg  and  Pc1  1n  the  gap  test  are  end  effects  of 
multidimensional  hydrodynamic  flow,  and  cannot  be  adequately  studied  within 
the  framework  of  simple  1-0  planar  theory.  However,  there  are  2-D  codes 
available  to  make  numerical  computations  Illustrating  the  flow  that  will  occur 
In  the  LSGT.  Computations  that  have  been  carried  out  on  the  donor/PMMA  system 
in  conjunction  with  particle  velocity  measurements  Indicate  that  ^he  lateral 
rarefaction  reaches  the  axis  of  the  PMMA  cylinder  (gap)  after  the  shock  has 
travelled  35  mm  (about  2/3  the  cylinder  diameter)  .  [The  arrival  of  the  release 
wave  Is  responsible  for  the  hump  which  appears  In  the  calibration  curve  of 
Fig.  18.  indeed  with  more  accurate  determination  of  P,  we  night  see,  Instead 
of  a  hump,  a  cusp  at  the  join  of  two  different  curves;  this  was  the  pattern 
found  In  the  u-x  curves.]  The  detailed  numerical  study  of  the  flow  caused  by 
the  shock  transmitted  to  the  acceptor  will  be  considered  after  an  experimental 
study  of  both  confined  and  unconfined  shocked  acceptors  has  been  described. 

53  - 

Kamegal,  M.  and  Erkman,  J.  0.,  "Numerical  Analysis  of  a  Diverging  Shock  Have 
In  Cylinders  of  Plexiglas,"  Fifth  Symposium  (International)  on  Detonation, 

ONR  ACR-184,  U.S.  Gov'nt  Print.  Office,  Washington,  DC,  1972;  pp  477-485. 
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Jaffe,  I.  and  Clalrmont,  A.  R. ,  Jr.,  "The  Effect  of  Configuration  and 
Confinement  on  Booster  Characteristics,"  NOLTR  65-33,  13  Apr  1965. 
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Cast  DINA  (diethanol  dinitrate)  was  the  acceptor  studied  for  Its  response 
to  shock55.  The  standard  LSGT  configuration  was  used  except  that  charge 
diameters  were  38.1  mm  (standard  for  unconfined  charges,  but  1.6  nin  larger 
than  the  standard  for  confined  charges)  and  152.4  mm  instead  of  139.7  mm  long. 

The  acceptors  were  Instrumented  with  continuous  resistance  wire  circuits  which 
measured  the  position  of  the  detonation  front  as  a  function  of  time.  Scopes 
were  triggered  by  a  pin  between  donor  pellets;  In  addition,  x-t  data  from  the 
LSGT  calibration  were  available.  The  data  were  reduced  to  obtain  run  length 
xs  and  delay  time  as  a  function  of  Initiating  pressure  P^.  Fig.  21  summarized 
the  data  xs  vs  for  the  confined  and  unconfined  charges,  and  shows  the  approach 
to  P  ^  for  this  experiment  In  each  case  (Pc1  value  Is  shown  by  vertical  lines 
at  the  asymptotic  values).  As  this  figure  shows,  the  curve  for  the  confined 
acceptor  conforms  exactly  to  that  expected  (Fig.  15).  The  curve  for  the 
unconfined  acceptor  is  Initially  (at  high  pressure  and  small  lengths  traveled 
by  the  shock)  coincident  with  the  curve  for  the  confined  charge.  This  merely 
means  that  confinement  has  no  effect  until  the  entering  curved  shock  front 
strikes  and  Interacts  with  the  confining  walls.  However,  the  two  curves  diverge 
rapidly  after  that,  and  a  cusp  appears  in  the  upper  one.  After  considering 
possible  causes  for  this  behavior  and  noting  a  similar  cusp  in  another  worker's 
unconfined  gap  testing  of  Comp  B,  we  concluded  that  it  is  an  effect  caused 
by  the  lateral  rarefaction  waves  and  perhaps,  In  small  part,  by  the  curvature 
of  the  entering  shock  front.* 

On  the  other  hand,  the  x  vs  Pj  curve  for  the  confined  (standard)  LSGT 

5  <  1 

shows  no  cusp  and  also  shows  a  linear  plot  x$  vs  P^  as  does  the  wedge  test 

which  will  be  described  later.  It  appears  that  the  measurements  were  completed 
before  the  arrival  of  lateral  rarefactions.  Hence  P„4  measured  in  the  standard 

FF -  - - — Cl 

33Pr1ce,  D.,  Toscano,  J.  P. ,  and  Jaffe,  I.,  "Development  of  the  Continuous  Wire 
Method.  Progress  Report  III,  Measurements  in  Cast  PINA,"  20  Apr  1967. 

♦It  Is  known  that  the  shape  of  the  entering  shock  wave  (planar  or  non-planar) 
affects  the  flow  and  presumeably  the  value  of  Ph  measured.  In  fact,  a  recent 
projectile  Impact  study  on  cast  Comp  showed  that  local  critical  energies 
for  Initiation  varied  with  the  area  impacted  and  were  generally  much  higher 
than  that  measured  (and  found  constant)  in  1-D  experiments.  That  suggested  a 
concept  of  critical  shocked  area  to  be  added  to  the  critical  parameters  of 
pressure  and  time. 

55Moulard,  H.,  "Critical  Conditions  for  Shock  Initiation  of  Detonation  by  Small 
Projectile  Impact,"  Preprints  of  the  Seventh  Symposium  (International)  on 
Detonation,  Vol.  1,  1981,  pp  182-189. 
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INITIATING  PRESSURE  (kbar) 


FIGURE  21  COMPARISON  OF  xs  VS.  Pj  FOR  CONFINED  AND  UNCONFINED  CHARGES. 
(REF.  55) 
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confined  LS6T  would  be  expected  to  approach  the  values  obtained  under  1-D 
planar  conditions;  there  would  be  no  similar  expectation  for  those  obtained 
with  unconfined  acceptors. 

Very  recently  hydrodynamic  calculations  have  been  made  of  the  flow  In  the 

57 

shocked  acceptor  of  the  LSGT  .  To  best  understand  the  results,  Information 
from  other  work  on  this  test  must  be  considered  at  the  same  time.  As  Fig.  17 
shows,  the  detonator  provides  essentially  a  point  Initiation  of  the  cylindrical 
pentolite  donor.  The  detonation  front  expands  spherically  as  it  travels  down 

CO 

the  cylinder.  Placesi  carried  out  a  hydrocalculation  (2D  Cyclone  code) 
of  such  a  divergent  flow  In  a  50/50  pentolite  cylinder  of  50.8  mm  dla;  his 
calculations  were  carried  out  to  an  a/d  of  1.7  (the  point  at  which  the  lateral 
rarefactions  eliminate  the  rear  rarefactions),  i.e.  well  beyond  the  *./d  of  1.0 
of  the  standard  donor.  Those  calculations  give  us  only  a  semi -quantitative 
picture  of  the  flow  in  the  standard  donor  because  the  doner  density  chosen  was 
1.65  g/err  Instead  of  the  standard  1.56  g/cm°  of  the  LSGT  donor.  In  particular, 
the  calculation  showed  that  the  detonation  front  in  the  cylinder  Is  that  of 
the  spherical  Taylor  expansion  wave  down  to  the  sonic  point.  As  Illustrated 
In  Fig.  22,  the  pressure  gradient  at  the  detonation  front  is  very  large.  Hence 
It  is  most  unlikely  that  the  Ideal  Pj  can  be  measured  or  that  a  corresponding 
transmitted  shock  Into  the  PMMA  can  be  attained.  This  fact  was  Ignored  In 
deriving  the  calibration  curve  of  Fig.  18  because  the  curve  was  forced  to 
terminate  (zero  gap)  at  the  pressure  that  would  have  been  transmitted  by  the 
one  dimensional  Pj. 

The  spherical  expansion  of  the  detonation  front  in  the  donor  resulted  In 
a  radius  of  curvature  (Rc)  equal  to  the  distance  traveled  by  the  front.  The 
shock  transmitted  to  the  PMMA  gap  shows  uninterrupted  continuation  of  the 
spherical  expansion  until  lateral  rarefactions  affect  it  at  about  Lg —  25  mm. 
There  Is  a  jump  In  the  pressure  at  the  interface;  the  pressure  transmitted  Into 
the  gap  Is,  of  course,  lower  than  the  donor  loading  pressure.  As  the  shock 

57 

Bowman,  A.  L. ,  Kershner,  J.  D.  ,  and  Mader,  C.  L.,  "A  Numerical  Model  of  the 
Gap  Test,"  LA-8408,  Oct  1980. 
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Placesi,  Jr.,  0.,  "Numerical  Hydrodynamic  Calculations  of  the  Flow  of  the 
Detonation  Products  from  a  Point-Initiated  Explosive  Cylinder,"  N0LTR  66-150, 
Jan  1967. 
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travels  from  the  PMMA  to  the  acceptor,  there  Is  another  pressure  jump  with  the 
amplitude  changed  again  according  to  the  Impedance  of  the  acceptor. 

Ref.  57  carried  out  calculations  with  the  2DE  (two  dimensional  reactive) 
code  In  which  the  Forest  Fire  burn  rate  was  used  for  all  explosives  except 
pentollte.  For  that  HE,  the  C-J  volume  burn  model  was  used.  With  this  model, 
the  divergent  detonation  wave  reached  a  steady  pressure  of  145  kbar  rather  than 
a  C-J  pressure  of  228  kbar  (Pj  of  215  kbar  was  used  for  Fig.  18),  and  the  error 
of  using  Pj  Instead  of  a  lesser  pressure  In  the  donor  loading  was  pointed  out. 

As  a  result,  the  computed  calibration  curve  of  Ref.  57  was  much  lower  at 

Lg  <  10  mm  (In  the  "nominal"  range  of  Fig.  18)  and  somev/hat  higher  at  Lg  >  10  mm 

than  the  LSGT  curve  of  Fig.  18.  Fig.  23  Illustrates  the  comparison. 

It  should  be  remarked  at  this  point  that  the  Taylor  wave  and  Cyclone  code 
calculations  of  Fig.  22  suggest  a  much  steeper  pressure  gradient  than  the 
approximately  constant  103-90  kbar  of  Fig.  23  at  Lg  <  10  mm.  Moreover,  experi¬ 
mental  data  for  both  shock  velocity  vs  distance  and  particle  velocity  vs 

49 

distance  curves  In  donor  shocked  PMMA  Indicate  steeper  slopes  .  In 

particular,  the  u-x  curve  obtained  by  the  electromagnetic  velocity  (EMV)  gage 

measurements  on  tetryl  loaded  PMMA  Included  10  experiments  at  x  <  10  mm. 

That  curve  Is  shown  In  Fig.  24  and  illustrates  the  steep  initial  slope.  (The 

Initial  point  was  for  plane  wave  boosting  of  the  tetryl  donor,  but  previous 

EMV  measurements  in  a  50.8  mm.dla  x  50.8  mm  long  tetryl  cylinder  had  resulted 

59 

In  the  same  values  for  both  plane  wave  and  detonator  point  initiation  . 

In  view  of  the  above,  It  seems  likely  that  both  curves  of  Fig.  23  are  In 
error  for  x  <  10  mm:  the  calibration  curve  of  Fig.  18  too  high  and  the 
computed  curve  of  Ref.  57  too  low  in  that  range.  Moreover,  the  EMV  data 
at  x  <  10  mm  for  pentollte  loading  of  the  PMMA  can  be  used  to  compute  the 
shock  amplitude  by  combining  Eq.  10. 

U  =  2.561  +  1.595  u 

with  Eq.  11  to  obtain 

P  *  1.185  (2.561  +  1.595  u)u 

^Edwards,  D.  J. ,  Erkman,  J.  0.,  and  Price,  D. ,  "The  Measurement  of  Particle 
Velocity  In  Pressed  Tetryl,"  NOLTR  72-83,  Aug  1972. 
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where  P  and  u  have  the  respective  units  of  GPa  and  mm/ys.  Then  if  the  data 
are  forced  to  terminate  on  the  experimental  point  at  x  *  0,  the  calibration 
curve  of  FI 9.  18  is  corrected  to  the  curve  derived  from  u-x  data  shown  In 
Fig.  25.  Both  from  shock  front  curvature  in  the  PMMA  and  inertial  effects  of 
the  A1  roll  gage,  errors  in  particle  velocity  measurements  would  be  expected 
to  result  In  values  below  rather  than  above  the  true  values.  Her.ce,  they  suggest 
that  the  Reference  57  curve  is  too  low  for  the  range  x  <  10  mm.  Moreover,  the 
EMV  measurements  result  in  P  values  that  fall  on  the  old  calibration  curve  for 
x  >  5  mm,  below  It  at  x  <  5  mm. 

As  the  Fig.  25  cm  ves  show,  the  Ref.  57  computation  and  the  experimental 
calibration  curve  agree  within  the  maximum  estimated  experimental  error  at 
x  >  6  mm.  Moreover,  whether  the  CO  volume  burn  of  the  pentolite  donor  represents 
the  exact  loading  of  the  PMMA  or  not,  the  computations  with  that  loading  provide 
much  information  about  the  shock  in  the  acceptor  of  the  LSGT.  In  our  standard 
test  with  steel  confinement,  it  shows  in  the  acceptor  an  entering  curved  shock 
front  which  rapidly  flattens  out  because  of  reinforcement  by  the  shock  inter¬ 
action  at  the  steel  walls.  It  also  shows  the  reaction  starting  at  the  walls 
and  not  at  the  charge  center  on  the  longitudinal  axis;  this  may  be  the  case 
or  may  be  an  artifact  of  the  sensitivity  of  the  Forest  Fire  code.  The  results 
for  the  unconfined  acceptor  are  those  expected:  increasing  curvature  of  the 
shock  front  with  increasing  penetration  of  the  charge,  and  reaction  starting 
on  axis.  Hence  with  the  pentolite  loading  of  Ref.  57,  the  computations  reinforce 
the  experimental  results  on  DINA:  the  standard  LSGT  (steel  confinement)  approxi¬ 
mates  the  1-D  condition  in  the  acceptor;  the  use  of  unconfined  acceptors  does 
not. 

Re*.  L8  was  devoted  to  the  LSGT  and  summarized  much  information  about  it. 
Some  of  that  information  was  repeated  (and  sometimes  updated)  above.  The 
remainder  will  be  summarized  very  briefly.  Increasing  temperature  decreases 
Pg;  the  effect  is  extreme  if  the  temperature  is  just  above  tc.  The  porosity 
effect  is  the  same  for  both  Group  1  and  2  HE:  Pg  increases  with  decreasing 
porosity;  the  increase  Is  particularly  sharp  as  d£  (and  dead-pressing)  Is 
approached  for  Group  2  materials.  See  Fig.  26  and  27.  The  particle  size 
effect  differs  for  granular  and  voidless  (cast  or  plastic  bonded)  charges. 

The  former  show  P  increasing  as  6  decreases;  the  latter  show  the  reverse 
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FIGURE  25  P  VS  x  CALIBRATIONS  FOR  S>'«ALL  GAP  VALUES 
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80/20 

AP/W 


FIGURE  27  LSGT  RESULTS  FOR  PRESSED  CHARGES  OF  AP  AND  AP/FUEL  MIXTURES. 

(GROUP  2).  (REF.  28).  (  O  AP  141,7 u;  □  AP  14B,8ju;  X  AP  146/Al  (7*1 1,96/6; 
A  AP  145/Al  (7*i),  90/10). 


%  TMD 

FIGURE  28  EFFECT  OF  ALUMINUM  ON  SHOCK  SENSITIVITY  OF  AP 
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trend.  Addition  of  wax  to  Group  1  explosives  increases  P^;  to  AP,  decreases 

It  (Fig.  27).  Addition  of  A1  to  cast  or  pressed  TNT  (at  95%  TMD)  Increases 

Pg;  A1  added  to  pressed  AP  appears  to  have  the  same  effect  (Increased  Pg) 

at  %  TMD  <  72,  but  the  opposite  one  at  80%  TMD  or  greater.  (See  Fig.  28). 

Although  we  have  concentrated  on  Pg,  we  have  mentioned  the  effect  of 

duration  several  times.  It  remains  to  consider  a  critical  energy  concept  and 

how  the  LSGT  values  conform  to  It.  The  existence  of  a  critical  initiation 

energy  has  been  proposed  many  times  in  the  past.  The  first  quantitative 

52 

proposal  was  made  by  Walker  and  Wasley  .  They  showed  that  for  several 
essentially  voidless  HE,  under  1-D  experimental  conditions,  there  is  a 
constant  energy/unit  area 

2 

E  *  c  Put  (12) 

f  o 

required  for  initiation.  Eq.  12  applies  for  a  square  pressure  pulse  of 

amplitude  P,  duration  t,  and  velocity  U  in  the  unreacted  explosive;  u  and 

p.  are  the  particle  velocity  and  density,  respectively.  Subsequently, 

0  60 

de  Longueville  et  al .  showed  that  the  concept  held  not  at  all  for  the 
liquids,  TNT  and  NM,  very  pocrly  for  granular  RDX  (pQ  =1.55  g/cm  ),  only 
approximately  for  plastic  bonded  HMX  and  RDX,  and  quite  well  for  PBX  9404 
and  cast  Comp  B3.  See  Fig.  29.  Despite  its  Inapplicability  to  liquid  and 
porous  HE  and  to  occasional  voidless  solid  HEs*,  the  energy  criterion  has 
served  as  a  useful  engineering  guide.  By  applying  it  to  essentially  voidless 
pressed  TNT1,  it  was  possible  to  obtain  a  quantitative  version  of  the  pressure¬ 
time  threshold  for  Initiation  outlined  in  Fig.  14.  Fig.  30  shows  the  result 

de  Longueville,  Y.,  Fauqulgnon,  C.,  and  Houlard,  H. ,  "Initiation  of  Several 
Condensed  Explosives  by  a  Given  Duration  Shock  Wave,"  Sixth  Symposium 
(International)  on  Detonation,  ONR  ACR-221,  U.S.  Gov.  Print.  Office, 
Washington,  DC,  1978;  pp  105-114. 

♦Plastic  bonded  TATB  behaves  like  the  liquids  with  regard  to  shock  initiation 
of  detonation®*. 

^Jackson,  R.  K. ,  Green,  L.  G. ,  Bartlett,  R.  H.,  Hofer,  W.  W.,  Kramer,  P.  E. , 
Lee,  R.  S.,  Nidick,  E.  0.,  Jr.,  Shaw,  L.  L. ,  and  Welngart,  R.  C.,  "Initiation 
and  Detonation  Characteristics  of  TATB,"  ibid  .  pp  755-765. 
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B  *  HMX-NYLON 
C  ■  RDX-TNT  (60/40) 

D  -  RDX-POLYBUTADIENE 
E  -  P8X  9404 
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FIGURE  29  ENERGY  THRESHOLD  VS  SHOCK  PRESSURE  (REF.  60) 
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for  pressed  TNT.  (The  same  procedure  indicated  that  cast  TNT  had  a  critical 
energy  about  four  times  as  large.)  As  Fig.  30  shows,  the  limit  curve  was 
drawn  from  the  critical  pressure  (taken  as  the  C-J  pressure  lasting  for  the 
reaction  time)  at  d  ■  dc,  and  the  value  measured  In  the  LS6T  falls  on  the 
curve  near  Its  lower  end.  Presumeably,  threshold  values  for  other  initiation 
tests  of  this  TNT  (e.g.,  gap,  wedges,  projectile,  and  booster  tests)  will 
also  fall  on  this  limit  curve.  In  other  words,  all  initiation  tests  values 
for  a  given  HE  should  be  related. 

The  Wedge  Test 

The  wedge  test  was  designed  as  a  plane-wave  system  so  arranged  as  to 
permit  continuous  optical  monitoring  of  the  wave  front  within  the  shocked 
explosive  .  A  typical  experimental  setup  is  shown  In  Fig.  31.  This  Is 
a  sophisticated  experiment  and  correspondingly  yields  much  Information 
from  the  streak  camera  record  synchronized  with  the  event.  The  emergence 
of  the  planar  shock  from  the  barrier  (a  material  of  known  Hugoniot),  Its 
tilt  (If  any),  and  the  position  of  Intersection  of  the  shock  waves  and  the 
Inclined  mylar  covered  face  of  the  wedge  are  recorded  as  a  function  of  time. 

The  Initial  pressure  of  the  shock  wave  entering  the  explosive  may  be  varied 
by  selecting  different  explosives  and  thicknesses  of  the  explosive  pad  and 
different  attenuators  and  thicknesses  of  the  barrier.  The  reduced  data  give 
P(x,t),  U(x,t),  run  distance,  and  delay  time. 

Although  most  gap  tests  use  a  cylindrical  charge  and  dc  is  defined  for  a 
cylindrical  charge,  there  Is  no  restriction  on  the  charge  geometry  used  In 
measuring  P^ .  The  choice  of  the  wedge  shape  (to  obtain  continuous  measurements) 
and  of  the  wedge  angle,  to  assure  that  release  waves  cannot  arrive  before  the 
shock  wave,  result  In  data  reproducing  that  expected  in  the  Interior  of  larger 
charges.  The  records  of  run  length  vs  initial  pressure  are  very  similar  to 
that  of  Fig.  21  for  the  standard  LSGT.  Jacobs  et  al.^6  concluded  that  their 
results  agree  reasonably  well  with  shock  initiation  work  on  cylinders  of  cross- 
sectional  area  comparable  to  the  face  area  of  the  test  wedge  provided  only 
axial  measurements  are  made  on  the  cylinders  and  the  initial  shocks  have 
comparable  pressure-time  profiles.  Fig.  21  Indicates  that  there  Is  more  apt 
to  be  agreement  between  wedge  and  the  standard  LSGT  ( 1 . e . ,  confined  charge) 
than  In  the  unconfined  case. 
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(SHOCK  ATTENUATOR) 


TEST 

WEDGE 


CAMERA 

SLIT 


SIDE  VIEW  FRONT  VIEW 

(TYPICAL  WEDGE ;  32  x  32  x  14  mm,  25°  WEDGE 
WITH  APEX  ANGLE  90°.) 


FIGURE  31  WEDGE-TEST  ARRANGEMENT.  (REF.  40) 


61 


A  comparison  of  wedge  test  results  with  the  NOL  LSGT  was  first  presented 
In  the  appendix  of  Ref.  50.  Subsequently  corrected  data  for  the  wedge  tests 
were  used  for  the  comparison  in  Ref.  28  together  with  better  Hugoniot  data 
available  at  that  time.  The  parameters  considered  in  reference  50  were: 
x$  Run  length,  distance  from  the  shocked  surface  to  the  plane  in 
which  steady-state  detonation  is  first  established. 
t$  Delay  time,  total  time  from  entry  of  shock  into  the  explosive 
until  the  time  at  which  steady-state  detonation  begins. 
a  Excess  delay  time,  amount  by  which  exceeds  the  time  required 
for  the  detonation  to  travel  distance  xg.  A  a  t$  -  xs/D. 

Pg  Pressure  (501  point)  in  PMMA  at  gap/acceptor  Interface. 

P.|  Shock  pressure. 

Pc1  Threshold  pressure  in  explosive  required  to  initiate  detonation 
in  LSGT. 

It  was  tentatively  concluded  that  a"1  vs  P,  t  1  vs  P,  and  x  -1  vs  P  were 
all  linear,  and  that,  in  the  last  case,  the  point  ( Pc ^ ,  0.02  mm  *)  fell  on 
the  curve.  Moreover,  it  was  pointed  out  that  log  x$  vs  log  P  was  also 
linear  and,  in  some  cases,  extrapolated  to  reasonable  values  of  P  for  x  ~ 

_  i  * 

reaction  zone  length  whereas  the  x,  vs  P  curves  did  not. 

5  50 

Because  the  averaae  set  of  data  consisted  of  only  three  points,  the 

46 

suggested  relationships  were  considered  tentative.  However,  Jacobs  et  al 
working  Independently  showed  linearity  of  xp  *  vs  P  for  some  nine  points  of 

5  fi? 

data  for  cast  Comp  B-3.  Somewhat  later  Ramsay  and  Popolato  published  linear 
log  x$  vs  log  P  curves  for  five  explosives.  Since  then  both  Los  Alamos 
Scientific  Laboratory  and  Lawrence  Livermore  Laboratory  have  used  this 
empirical  relationship  for  evaluating  shock  sensitivities. 

Data  for  confined  DINA  (LSGT-Fig.  21)  gave  a  linear  curve  xs_1  vs  P 
whereas  data  for  unconfined  DINA  showed  a  sharp  change  in  slope  and  non¬ 
linearity  over  the  low  pressure  portion  of  the  xs"*  vs  P  curve. 

Ramsay,  J.  B.  and  Popolato,  A.,  "Analysis  of  Shock  Wave  and  Initiation  Data 
for  Solid  Explosives,"  Fourth  Symp.  (Int.)  on  Detonation,  0NR,  ACR-126, 

U.S.  Gov.  Print.  Office,  Washington,  D.C.,  1967,  pp  233-238. 
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Data  for  the  confined  DINA  were  also  examined  In  the  other  empirical 

relationships  that  have  been  mentioned  (log  x_.  vs  log  P,  log  x_  vs  log  A, 

•  1  1  '  * 
t$  vs  P,  i  vs  P).  All  showed  some  curvature  when  all  eight  pairs  of  data 

were  Included.  However,  If  the  first  two  pairs  (l.e.,  at  P  values  of  42.8 

and  49.5  kbar)  were  discarded,  the  remaining  six  showed  linearity  for  all 

the  other  relationships  mentioned.  Thus  for  the  confined  DINA  and  presumably 

any  comparable  data  obtained  In  the  LSGT  set-up,  the  most  successful  linear 

relation  Is  that  of  x  vs  P. 

-1 

This  same  relation  (x$  vs  P)  has  been  used  for  reexamination  of  the 
wedge  test  data.  The  transition  data  for  cast  pentollte  and  DINA  are  plotted 
In  Figure  32;  those  for  the  remaining  cast  materials,  in  Figure  33.  The  data 
for  pentollte  and  DINA  have  been  plotted  separately  (1)  to  emphasize  the 
similarity  between  the  wedge  test  results  and  those  from  the  regular  LSGT 
configuration  on  similar  explosives  and  (2)  because  these  curves  cross  those 
of  Figure  33  for  the  other  TNT  based  castings.  Both  the  pentollte  and  DINA 
data  (Xj-1  vs  P)  seem  to  terminate  at  P^.  The  apparent  crossing  of  these 
curves  near  P^  Is  probably  fictitious.  The  P^  values  of  11.9  (pentollte) 
and  11.0  kbars  (DINA)  are  not  significantly  different  and  could  be  reversed 
by  small  errors  In  the  measurements  or  the  Hugonlots  used  or  both.  Moreover, 
unconfined  pentollte  showed  Pci  of  12  -  14  kbar  as  compared  to  15  kbar  for 
unconfined  DINA,  l.e.,  the  reverse  of  the  relative  sensitivities  of  the 
confined  charges.  If  the  run  length  for  the  same  shock  amplitude  Is 
considered  a  measure  of  the  shock  sensitivity,  then  cast  pentollte  is  more 
shock  sensitive  than  cast  DINA,  probably  over  the  entire  pressure  range. 
However,  both  curves  of  Figure  32  intersect  some  of  the  curves  of  Figure  33. 
Hence  on  the  run-length  criterion,  reversals  in  relative  sensitivity  will 
occur  between  the  higher  and  lower  pressure  ranges,  e.g.,  for  pentollte  and 
Comp  B-3  and  for  DIMA  and  Octol  65/35. 

In  Figure  33,  all  five  cast  explosives  show  a  satisfactory  linear 
extrapolation  of  xs_1  vs  P  to  P  ^ ,  save  possibly  Comp  P.  Of  course  the 
data  for  Comp  B-3  Include  only  the  higher  pressure,  shorter  run  length 
data.  Reference  44  gives,  in  uddition,  lower  pressure,  longer  run  length 
data  that  resulted  In  an  extrapolation  to  28  kbar  instead  of  the  present 
22  kbar.  This  difference  is  the  order  of  magnitude  of  errors  In  P  ^ 
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N  LENGTH 


SHOCK  PRESSURE,  P  (kbor) 


FIGURE  32  SHOCK  TO  DETONATION  TRANSITION  FOR  CAST  PENTONITE  AND  DINA. 

(O  PENTONITE,  WEDGE  DATA;  □  DINA,  REGULAR  LSGT  CONFIGURATION; 
XPcj  VALUE).  (REF.  28) 


FIGURE  33  SHOCK  TO  DETONATION  TRANSITION  FOR  OTHER  TNT  BASED  CAST  EXPLOSIVES. 

(O  COMP  B-3;  □  OCTOL  65/35;  A  COMP  B;  0  CYCLOTOL  ?o/25;  •  TNT;  X  Pci  VALUE).  (REF.  28) 
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and  possibly  of  the  lower  shock  pressure  measurements.  It  might  also  be 
caused  by  an  unmonitored  duration  effect  In  the  action  of  the  shock  wave; 
such  an  effect  becomes  more  important  at  the  lower  amplitudes.  Whatever 
the  cause (s),  the  discrepancy  cannot  be  resolved  without  further 
experimental  work. 

By  the  run-length  criterion,  the  shock  sensitivities  of  Figure  33,  in 
descending  order,  are  Comp  B-3  (RDX/TNT,  60/40),  Octol  65/35,  Comp  B 
(RDX/TNT/Wax,  60/40/1),  Cyclotol  75/25,  and  TNT.  This  ordering  of  cast 
explosives  can  be  reconciled  with  composition  effects  if  the  role  of 

cq 

particle  size  in  affecting  shock  sensitivity  is  recognized.  Seely 
demonstrated  in  1963  that  shock  sensitivity  increased  with  increasing 
particle  size  If  air  was  the  continuous  medium,  but  increased  with  decreasing 
particle  size  if  the  continuous  medium  was  a  condensed  one.  Cast  explosives 
fall  in  the  latter  category  and  should  show  Increasing  sensitivity  with 
decreasing  particle  size.  RDX  Class  F  is  recommended  for  preparing  Comp  B-3; 
Class  4  for  Comp  B;  and  Class  D  for  cyclotol s.  The  relative  weight  mean 
diameters  are  58,  130,  and  740u  respectively  for  Classes  F,  A,  and  D.  It  Is 
evident  that  the  difference  in  the  RDX  particle  size  can  account  for  the  fact 
that  cast  Comp  B-3  is  more  shock  sensitive  than  cast  Comp  B.  The  difference 
between  Classes  A  anti  D  RDX  can  also  account  for  the  cast  cyclotol  ( 75%  RDX) 
being  less  sensitive  than  the  cast  Comp  B  (60%  RDX).* 

Finally,  Figures  32  and  33  reemphasize  the  inadequacy  of  an  attempt  to 
characterize  shock  sensitivity  by  measurement  at  only  the  50%  point.  The 
curves  xs_1  vs  Pi  not  only  demonstrate  differences  in  sensitivity  much 
better  than  LSGT  values  alone  can  do;  they  also  show  reversals  in  behavior 
that  cannot  be  predicted  from  the  LSGT  values  unless  an  instrumented  LSGT 

go 

Seely,  L.  B. ,  "A  Proposed  Mechanism  for  Shock  Initiation  of  Low  Density 
Granular  Explosives,"  Proc.  Fourth  Elec.  Initiation  Symposium,  Franklin 
Institute,  Phila.  1963;  Paper  27  of  Rept.  EIS-A-2357. 

*When  pressed  cyclotol  75/25  and  pressed  Comp  B  are  compared  at  96.4%  TMD  In 
the  LSGT  their  respective  50%  values  of  Pg  are  12  and  18  kbar.  By  going  to 
highly  pressed  charges,  the  Initial  particle  size  effect  has  been  reversed 
and  the  chemical  effect  (75  vs  60%  RDX)  is  still  present.  Hence  the  relative 
sensitivity  has  not  only  been  reversed  (compared  to  ratings  of  the  cast 
materials),  but  the  difference  is  significant. 


65 


NSWC  TR  80-339 


(such  as  that  of  Fig.  32)  Is  used.  Even  Figures  32  and  33  are  incomplete 

descriptions  of  shock  sensitivity.  As  discussed  earlier,  a  limit  curve  in 

the  pressure-time  plane  exists  for  the  shock  initiation  of  detonation  In 

each  charge,  and  a  similar  limit  curve  probably  exists  for  the  shock 

initiation  of  burning  and  sub-detonation  reactions.  Liddiard^  has  shown 

that  threshold  conditions  for  initiating  such  reactions  can  be  defined  and 

measured.  Moreover  such  shock  induced  reactions  are  also  an  important  aspect 

of  a  material's  sensitivity  to  shock. 

Relations  Between  Results  of  Various  Gap  Tests 

In  accord  with  the  idea  that  all  P  .  values  measured  in  different  tests 

ci 

should  lie  on  curves  such  as  those  of  Figs.  14,  16,  30  (for  a  specific  explosive) 

and  should  be  related  to  each  other,  correlations  between  different  gap  tests 

were  investigated.  It  has  been  shown  that  the  MOL  small  scale  gap  test 

(SSGT)  results  correlate  well  with  the  LSGT  results.  Surprisingly,  it  has 

been  shown  'hat  LSGT  results  for  confined  and  unconfined  cast  HE  acceptors 

also  correlate.*  Hence  a  relationship  might  be  expected  between  any  other 

well  standardized  shock  initiation  tests  carried  out  on  identical  explosive 

samples.  Such  comparisons  with  results  of  other  laboratories  have  to  be  made 

with  502  point  gap  length  (L  )  values  because  of  uncalibrated  tests  or  differences 

y 

in  calibration  data  used.  In  contrast  to  the  use  of  Pg  rather  than  Pc- ,  the 
use  of  gap  length  distorts  the  true  sensitivity  scale,  particularly  in  the 
higher  sensitivity  (lower  P^ )  range.  For  two  tests  of  quite  similar  dimensions, 
one  would  hope  both  would  be  affected  in  the  same  manner;  if  so,  any  distortion 
would  not  conceal  the  relation  between  the  tests  when  comparing  502  gap  values. 

LASL  has  a  LSGT,  unconfined,  of  41  mm  di a  ;  this  compares  to  the  NOL 
LSGT,  unconfined,  of  38  mm  dia  which  correlates  wel 1  with  the  standardized 
LSGT  values.*  Hence  despite  different  donors  and  different  attenuators  some 
correlation  between  test  results  on  the  same  HEs  would  be  expected. 

p - 

Liddiard,  T.  P. ,  "The  Initiation  of  Burning  in  High  Explosives  by  Shock 
Waves,"  Fourth  Symposium  (International)  on  Detonation,  ONR  ACR-126, 

U.S.  Gov.  Print.  Office,  Washington,  D.  C. ,  1967;  pp  487-495. 

*lhe  correlation  holds  only  for  the  502  point  gap  length  Lg ;  if  P?  is  used,  the 
data  show  a  jump  at  x  =  35  mm,  the  point  at  which  rarefactions  reach  the  axis 
of  the  PMMA  cylinder. 

„  -  '  \ 
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Table  2  contains  the  reported  data  for  seven  cast  explosives  tested  at 
both  LASL  and  NSWC.  All  LASL  data  came  from  Reference  23;  NSWC  data  are  from 
Reference  28  (Tests  made  after  8  Mar  1974  are  designated  by  test  number  in 
the  table.)  Differences  in  composition,  density,  and,  where  known,  preparation 
are  shown.  LASL  In  testing  eight  batches  of  ca  cyclotol  had  a  range  of 
165-180  (average  value  of  172  ±  8)  inches  x  102;  NSWC  recorded  tests  for  five 
batches  of  cast  50/50  pentolite  with  a  range  of  272-284  (average  value  of 
278  ±  6)  in.  x  10^.  The  two  laboratories  seem  comparable  in  reproducing  their 
values  from  batch  to  batch  of  cast  (TNT  matrix)  HE.  Figure  34  shows  a  plot 
of  the  LSGT  data.  The  two  sets  are  definitely  related  with  L  (NSWC)  consistently 
greater  than  Lg(LASL).  The  exact  nature  of  the  relationship  can  only  be 
established  with  a  greater  number  of  explosives  of  the  same  composition  and 
preparation  tested  by  both  laboratories.  For  example,  the  data  of  Figure  34 
might  be  better  fit  with  a  curved  than  a  straight  line;  the  differences  in 
composition  and  procedures  (Table  2)  will  certainly  affect  the  final  result. 

Even  with  the  present  data,  there  are  no  reversals  ir  sensitivity  ratings. 

(There  is  an  apparent  reversal  between  Comp  B-3  and  octol ,  but  these  two  HE 
have  the  same  LSGT  value,  within  experimental  error,  in  both  sets  of  data. 

Hence  the  apparent  reversal  is  not  considered  significant.) 

The  situation  for  pressed  explosives  is  more  complex.  Table  3  contains 
the  data  which  arc  plotted  in  Figure  35.  With  three  exceptions,  noted  in 
both  table  and  gr„,  i,  all  the  explosives  were  pressed  at  25°C.  In  general, 
explosives  pressed  at  different  temperatures  would  not  be  comparable.  Figure 
35  seems  to  suggest  that  there  is  a  relationship  between  the  two  sets  of  test 
values  for  si::  of  eight  pressed  HE.  The  other  two  HE  are  DATES  and  Comp  A-3. 

The  difference  of  the  former  might  be  attributed  to  the  comparison  of  a  cold 
pressed  charge  with  a  hot  pressed  charge.  In  the  case  of  Comp  A-3,  the  NSWC 
charges  were  prepared  from  plant  production  HE  which  we  have  found  in  earlier 
work  gave  rather  erratic  LSGT  results.  That  may  also  have  been  the  case  at 
LASL;  as  noted  in  the  table,  higher  results  were  reported  in  Reference  23 
for  other  A-3  charges  of  nearly  the  same  density. 

Although  the  large  departures  of  DATB  and  Comp  A-3  from  the  average 
behavior  may  be  explained  as  above,  the  different  order  of  sensitivity  ratings 
cannot  be  ignored.  The  ratings  are  shown  in  the  lower  part  of  Table  3.  It  is 
possible  that  hot  pressing  of  DATB  so  sensitizes  it  that  it  seems  equivalent 
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100  200  300 

LASL  GAP  (IN.  x  102) 

FIGURE  34  COMPARISON  OF  LSGT  VALUES  FOR  CAST  HE 
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FIGURE  36  COMPARISON  OF  LSGT  VALUES  FOR  PRESSED  HE  (>90%  TMD) 
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to  Explosive  D.  It  Is  also  possible  that  an  atypical  value  for  Comp  A-3  Is 
responsible  for  the  Inversion  of  the  TNT/A3  values  In  the  LASL  data  compared 
to  those  of  NSWC.  The  large  difference  between  hot  and  cold  pressed  TNT, 
seen  In  the  LASL  data,  does  not  appear  In  the  NSWC  data.  That  might  well  be 
because  little  of  the  heat  from  the  water  medium  of  the  Isostatic  press  was 
transmitted  thru  the  rubber  meld  to  the  TNT.  In  addition  to  all  these 
possibilities,  Initial  particle  size  also  has  an  effect  on  LSGT  values. 

At  the  present  time,  It  does  not  seem  warranted  to  assume  a  reliable 
relation  between  the  two  sets  of  test  values  for  pressed  HE.  This  view 
Is  reinforced  by  the  definite  lack  of  correlation  In  the  case  of  porous  TNT 
(pressed  at  25°C).  The  data  are  listed  In  Table  4  and  plotted  In  Figure  36. 
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FIGURE  36  COMPARISON  LSGT  VALUES  FOR  POROUS  TNT 
(CHARGES  PRESSED  AT  25°C) 
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INTERDEPENDENCE  OF  dc  AND  P  ^ 

From  the  information  now  available  on  critical  diameter  for  detonation 
and  critical  pressure  for  initiation  it  seems  evident  that  the  value  of  either 
of  these  quantities  Is  one  of  the  factors  determining  the  other.  However,  It 
Is  not  the  dominating  factor;  hence  no  general  relationship  between  only  dc 
and  P  j  exists. 

Fig.  37  Illustrates  the  qualitative  effects  on  dc  and  P  ,  of  Increasing 
the  value  of  a  number  of  variables  defining  the  experimental  conditions.  As 
It  shows,  confinement,  temperature,  and  dilution  with  wax  cause  parallel 
effects.  Porosity  causes  similar  effects  In  Group  2  HE  (e.g.,  AP)  but  opposite 
effects  In  Group  1  HE  (e.g.,  TNT).  Particle  size  causes  similar  effects  In 
voidless  HE,  but  opposite  ones  in  porous  HE.  Addition  of  A1  also  seems  to 
cause  opposite  effects  over  most  of  the  range.  Out  of  11  pair  comparisons, 
there  were  four  reversals;  this  fact  alone  rules  out  a  general  relationship 
between  only  these  two  critical  parameters. 

On  the  other  hand,  a  continuous  variation  of  dc  with  exists  for  a 
single  material  in  various  forms.  Fig.  38  for  TNT  is  a  good  illustration 
of  this  because  this  material  is  castable.  From  very  porous  granular  charges 
through  hot  pressed  charges  to  various  forms  of  castings  offers  a  large  variety 
of  physical  states  as  well  as  introducing  the  variable  of  physical  heterogeneity 
not  previously  mentioned.  Although  this  quantity  cannot  be  measured  precisely, 
we  know  that  it  is  greatest  at  the  left  of  Fig.  38  and  least  at  the  right. 
Moreover,  we  see  the  trend  of  the  curve  dc  vs  P  ^  reverse  itself  as  the  porous 
charges  are  replaced  by  voidless  ones.  Thereafter  the  critical  values  vary 
in  the  same  way:  increasing  with  decreasing  heterogeneity  (increasing 
homogeneity).  The  effect  of  degree  of  heterogeneity  is  attributed  to  its 
effect  on  hot  spot  formation. 

Fig.  39  shows  the  relations  dc  vs  Pg  for  various  granular  Group  1  HE  - 
all  pure  explosives.  The  dashed  lines  connect  all  the  values  at  Bh%  TMD 
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FIGURE  37  EFFECT  OF  VARIABLE  CHANGE  ON  VALUES  OF  dc  AND  Pg 
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to  demonstrate  that  there  is  no  smooth  curve  (at  constant  porosity)  for  all 
explosives  although  the  curves  for  single  explosives  are  certainly  smooth. 

The  dashed  line  also  illustrates  that  the  physical  form  of  NQ  (low  bulk  density, 
and  two  different  lots  of  high  bulk  density  material)  had  very  little  effect 
on  the  Pg  value  at  85%  TMD  compared  to  a  relatively  large  effect  on  dc< 

Fig.  40  shows  comparable  d£  vs  Pg  relations  for  various  granular  Group  2 
HE,  AP  and  some  of  Its  mixtures.  In  this  case,  the  80?  TMD  points  are  marked 
to  demonstrate  again  that  there  Is  no  single  curve  encompassing  all  the  HE, 
despite  the  fact  that  each  one  has  a  smooth  curve.  Here  the  particle  size  of 
AP  at  ~  55%  TMD  (Initial  points)  has  more  effect  on  dc  than  on  Pg,  but  at  80% 

TMD  the  relative  effects  are  reversed. 

So  far,  the  relations  considered  have  been  for  porous  HE.  Ref.  28  contains 
some  data  for  cast  and  highly  pressed  explosives,  i.e.  essentially  voidless 
HE.  No  correlation  between  P  ^  and  dc  was  found  for  these  materials  although 
it  is  general  knowledge  that  systematic  variation  of  component  concentration 
in  multicomponent,  voidless  HE  frequently  results  in  such  a  correlation,  but 
only  for  each  series;  it  is  therefore  much  like  the  single  curves  of  Figs.  39 
and  40. 

Recently  it  was  suggested  that  a  linear  relationship  exists  between  In  dr 

24  c 

and  Lg  for  voidless  charges  .  In  view  of  the  relations  found  earlier  between 

the  NSWC  and  LASL  LSGT  values,  the  suggestion  was  investigated  for  each  set 

of  values  separately.  The  LASL  and  NSWC  data  are  given  in  Tables  5  and  6, 

respectively. 

The  d  values  of  Table  5  are  for  the  most  part,  LASL  values.  One  exception 

•  c  36 

is  that  for  TNT.  As  discussed  in  Ref.  21,  Cybulski  et.  al .  chose  unfortunate 

names  for  the  various  TNT  castings  they  prepared  and  tested.  To  quote  them: 

"(1)  Poured-clear. . . 

"(2)  Poured-cloudy.  The  molten  TNT  is  allowed  to  cool  to  the  region  of 
its  set-point  with  continuous  stirring  until  a  slurry  is  formed  which  Is  then 
poured  into  molds. 

"(3)  Creamed.  The  cooling  melt  is  agitated  with  a  mixer  of  the 
perforated  disk  type;  when  the  slurry  has  formed,  10%  by  weight  of  powdered 
TNT  of  such  fineness  that  all  passed  through  a  no.  52  BSS  sieve  (300  y 
opening)  is  added  slowly;  after  further  agitation,  the  thickened  melt  Is 
poured," 
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Procedure  No.  2  is  that  which  NSWC  uses  to  produce  "created  cast  TNT", 
its  customary  cast  TNT.  Consequently,  in  Table  4  Cybulski's  value  for 
"poured-cloudy"  TNT  Is  used  instead  of  the  value  of  his  so-called  "creamed" 

TNT.  Unfortunately,  Campbell  &  Engel ke14,  after  treating  the  data  for  the 
latter,  labeled  it  "creamed"  without  any  explanation.  It  Is  not  surprising 
that  the  value  of  dc  they  obtained  was  about  the  same  as  that  they  derived 
from  the  Stesik  and  Akimova  data  obtained  from  TNT  cast  In  water-cooled  molds 
(hence,  rapidly  cooled  to  produce  very  small  grains);  it  is  also  about  a  third 
of  the  values  for  creamed  cast  TNT  quoted  in  Table  4. 

In  Figure  41,  the  data  of  Table  5  are  plotted  In  d  vs  LSGT  value  L  . 

V  M 

Squares  are  used  for  pressed  charges,  of  which  there  were  three,  and  circles 
for  cast  charges.  It  is  evident  that  no  single  line  fits  all  the  data.  If 
the  points  for  tritonal  and  X0219  are  omitted,  a  straight  line  fits  the  six 
remaining  data  points  moderately  well.  However,  there  is  no  justification 
for  omitting  two  points. 

When  all  NSWC  LSGT  values  are  used  for  a  similar  comparison,  the  results 
are  somewhat  different  --  possibly  because  the  data  arc  more  extensive.  Table  6 
contains  the  available  data  which  are  plotted  in  Fig.  42.  Interestingly  enough, 
a  single  line  would  encompass  (with  moderate  scatter)  all  six  pressed  charges 
as  well  as  two  plastic  bonded  castings.  However,  all  of  the  TNT  based  castings 
would  require  another  quite  distant  parallel  line  and  even  then  have  several 
points  off  by  an  appreciable  arrounl .  Hence,  neither  LASL  nor  NSWC  data  indicate 
a  general  relationship  between  In  dc  and  Lg. 

For  some  time  now  a  group  at  Ballistics  Research  Laboratories,  studying 
Initiation  of  explosives  by  projectile  impact,  have  treated  ignition  (first 
flash  of  light  observed)  and  buildup  of  reaction  to  propagation  (monitored 
by  pressure-time  and  velocity-time  data  at  the  fronts)  as  separable  phenomena^9"7* 

6  Frey,  R. ,  Howe,  P.  ,  Timble,  J.,  and  Melani,  G.,  "Initiation  of  Explosive 
Charges  by  Projectile  Impact,"  ARBRL-TR-02176,  Oun  1979. 

70Taylor,  B.  C.  and  Ervin,  L.  H.,  "Separation  of  Ignition  and  Buildup  to 
Detonation  in  Pressed  TNT,"  Sixth  Symposium  (International)  on  Detonation, 

ONR  ACR-221,  Gov.  Print.  Office,  Washington,  D.C.,  1978,  pp  3-10. 

71Howe,  P.,  Frey,  R. ,  Taylor,  B.,  and  Boule,  V.,  "Shock  Initiation  and  the 
Critical  Energy  Concept,"  ibid,  pp  11-19. 
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They  have  shown  that  some  variables  have  opposite  effects  on  the  two,  e.g.  a 

large  particle  size  Increases  the  ease  of  Ignition  but  decreases  the  ability 

72 

of  the  reaction  to  build  up.  E.  Lee  et  al.  have  treated  Initiation  of 

heterogeneous  explosives  theoretically  as  an  ignition  followed  by  growth  of 

combustion  at  the  Ignlteu  "hot  spots".  They  have  had  success  with  PBX  9404, 

a  voldless  plastic  bonded  HE,  In  quantitative  agreement  between  predicted  and 

73 

experimental  values  of  P-t,  u-t  In  the  growth  region  . 

If  the  gap  test  Is  carried  out  at  d  a  3d  ,  It  will  be  under  experimental 

0 

conditions  so  favorable  to  reaction  growth  and  propagation  that  It  can  be 
considered  a  test  of  Ignltablllty  under  shock.  Similarly,  If  critical  diameter 
measurements  are  made  with  a  donor  considerably  more  energetic  than  the  test 
explosive  to  assure  Ignition  under  shock,  the  determination  of  d  can  be 
considered  a  qualitative  Inverse  measure  of  the  ease  of  reaction  build  up  and 
propagation.  The  above  is  an  oversimplified  picture,  of  course.  If  Ignition 
Is  by  hot  spot  formation,  as  Is  currently  believed,  then  the  more  hot  spots 
formed,  the  stronger  the  Ignition  and  the  buildup  and  propagation.  But  that 
assumes  the  usual  transport  processes  operating  throughout.  Obviously,  If 
detonation  Is  achieved,  the  propagation  will  be  hydrodynamic  and  Independent 
of  the  Initial  density  of  preformed  hot  spots.  Thus  hydrodynamic  propagation 
In  Group  1  HE  vs  high  pressure  burning  surface  reactions  In  Group  2  HE  could 
explain  the  porosity  reversal  of  trends  of  df  vs  In  Fig.  37.  Hot  spot 
formation  In  porous  HE  could  explain  the  reversal  of  the  particle  size  effect 
In  Fig.  37  as  well  as  that  caused  by  addition  of  Al . 


- - 

Green,  L.,  Nldlck,  E.,  Lee,  E. ,  and  Tarver,  C.,  "Reactions  in  PBX-9404 
from  Low  Amplitude  Shocks,  "  Behavior  of  Dense  Media  Under  High  Dynamic 
Pressures,  CEA  Paris,  1978;  pp  115-126.. 

73 

Lee,  E.,  "Initiation  of  Explosives  and  Propellants,"  presentation  at 
NSWC/WOL,  5  Jun  1980. 
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SUMMARY 

The  value  of  the  critical  diameter  d  is  determined  by  charge  temperature, 

c 

initial  particle  size,  density,  confinement,  composition,  and  the  nature  of 
the  HE  --  particularly  its  C-J  and  von  Neumann  pressures.  Sinte  detonation 
must  be  achieved  in  measuring  d  ,  a  shock  must  be  formed  and  the  snock 

v 

sensitivity  of  the  material  must  have  a  role. 

The  value  of  the  critical  Initiating  pressure  P  ^ ,  an  inverse  measure 
of  shock  sensitivity,  is  a  function  of  charge  temperature,  initial  particle 
size,  density,  test  dimensions,  confinement,  composition  and  the  nature  of 
the  material  --  particularly  its  d  and  function  d/d  . 

The  "shock  sensitivity"  of  an  explosive  cannot  be  described  by  a  single 
value  such  as  P^.  Shock  sensitivity  encompasses  all  of  the  complex  reactions 
of  the  explosive  to  many  different  shock  conditions.  Such  reactions  are 
manifested  by  the  phenomenon  of  failure  diameter  for  propagation  of  detonation 
as  well  as  by  the  critical  Initiation  pressures  for  shock  to  detonation 
transition  observed  in  the  numerous  gap,  booster,  projectile,  and  wedge  experi¬ 
ments.  A  complete  description  of  an  explosive's  sensitivity  to  shock  should 
encompass,  therefore,  not  only  the  thresholds  for  initiation  under  transient 
conditions  but  also  the  threshold  for  propagation  under  steady  state 
conditions.  Moreover,  it  should  account  for  the  large  effect  of  the  physical 
state  of  the  explosive  on  both  thresholds.  The  relationships  between  these 
many  aspects  of  shock  sensitivity  are  complex  and  cannot  be  simply  stated. 

When  P  ^  is  measured  in  two  different  gap  tests,  the  twc  sets  of  values 
obtained  from  a  group  of  the  same  HE  charges,  identically  prepared,  correlate 
well  provided  geometry  and  conditions  of  the  two  tests  are  similar.  Wedge  test 
results  and  shock-to-detonation  transition  studies  on  cylindrical  charges  of 
cast  DINA  showed  strong  similarities.  It  appears  that  the  standard  LSGT  values 
lie  on  the  low  pressure  end  of  the  curve,  reciprocal  run  length  vs  initiating 
pressure,  defined  by  the  wedge  test  data. 
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Although  dc  and  Pcj  are  obviously  Interdependent,  this  fact  does  not 
dominate  the  effects  of  other  variables  to  produce  a  general  relation  between 
only  these  two  parameters.  Instead  of  d  ,  It  Is  sometimes  more  convenient  to 

V 

use  critical  density  at  a  given  diameter  for  porous  explosives.  The  emphasis 
on  dc  and  Pc,1  is  such  that  other  critical  parameters  (e.g.  critical  density, 
particle  size  and  temperature)  are  apt  to  be  Ignored. 
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